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SYMBOLS 


1 9 to R 

overheating parameter, — ^ 

R - R 
"w r 

a overheating parameter, 

W x\ 

r 

c_p skin friction coefficient 

specific heat at constant pressure 

d diameter of hot wire 

E bridge voltage 

f frequency 

H heat loss from wire 

h enthalpy 

I current 

d £n R 

K wire parameter, - — - — — 

* 9 d In T 

w 

k thermal conductivity 

£ wire length 

length of separated flow 
M Mach number 

m boundary-layer mass flow 

Nu Nusselt number 

n x direction cosine of sound wave 

p pressure 

p t P pitot pressure 

total he at- transfer rate 
molecular heat-transfer rate 

R g resistance in series with wire (includes cable resistance) 

R resistance of wire 

w 



Re Reynolds number 

Re^ wire Reynolds number, pud/y^ 

r coordinate defined in figure 2 

r fluctuating diagram parameter 

5 fluctuating diagram parameter 

S* fluctuating diagram parameter 

T temperature 

T, . total stress tensor 

ij 

t time 

u velocity component defined in figure 2 

-v ->• ->■ 

V total velocity, ui + vj + wk 

v velocity component defined in figure 2 

W power to wire 

w velocity component defined in figure 2 

x coordinate defined in figure 2 

y coordinate defined in figure 2 

a 1/(1 + [( Y -1)/2]M 2 ) 

6 (y-l)M 2 a 

y ratio of specific heats (l.UO for present study) 

6 “boundary-layer thickness 

6* boundary-layer displacement thickness 

Ae^ fluctuation sensitivity 

n recovery factor, T r / T t 

tt sound pressure mode amplitude, 100 p f /YP 

0 overheating parameter, T^/T^; or coordinate defined in figure 2 

p density 

X wavelength 

a vorticity mode amplitude, 100 u’/u 


vi 


t entropy mode amplitude, 100 T f /T 

r molecular stress tensor 

ip — T 

w r 

overheating parameter, — 

<() yaw angle of wire relative to local flow angle 

y absolute viscosity 

\il second coefficient of viscosity 

a) local flow angle relative to nozzle center line 

wi flow turning angle across incident shock wave 

Subscripts 

c convected in direction of flow 

e edge of viscous effects 

r recovery conditions 

t local total conditions 

wall evaluated at the wall 

w pertaining to wire 

00 conditions ahead of incident shock wave 

Superscripts and Other Symbols 
( ) f fluctuating quantity 

(") time-averaged quantity 

<()> rms of fluctuating quantity 



THE BEHAVIOR OF A COMPRESSIBLE TURBULENT BOUNDARY LAYER IN A 

t 

SHOCK-WAVE-INDUCED ADVERSE PRESSURE GRADIENT 
William C. Rose 
Ames Research Center 

SUMMARY 

A review of previous analytical and experimental investigations of the 
behavior of a compressible, turbulent boundary layer in a shock-wave-induced 
adverse pressure gradient is presented. 'This review led to the following 
assessment of the present status of research on such boundary- layer flows: 

!• Present analytical methods for describing flow in such pressure 
gradients are incapable of properly accounting for the turbulent mixing prop- 
erties in the boundary layer. 

2. These mixing properties change throughout the flow and they signifi- 
cantly affect the boundary- layer flow downstream of the shock-wave - 
boundary-layer interaction. 

3. Data on which to base improved models of the turbulent mixing charac- 
teristics for such compressible flows are not available. 

The purpose of the present investigation was to obtain mean- and 
fluctuating- flow data, including the complete Reynolds stress tensor and tur- 
bulent mass- and heat-transfer rates, upstream of, within, and downstream of a 
severe adverse pressure gradient induced by an externally generated shock wave. 
This body of data can be used to assess the validity of analytical prediction 
methods, and to aid in formulating improved turbulent mixing models. 

The differential conservation equations of mass, momentum, and thermal 
energy are written in terms of mean and fluctuating quantities and then time 
averaged. These equations yield the turbulent transport terms that affect the 



time-averaged (mean) flow. The present work reports the results of an experi- 
mental investigation to determine each of these turbulent transport terms. 

The experimental investigation was conducted in an axially symmetric flow 
at a nominal free-stream Mach number of h and a Reynolds number based on 
boundary- layer thickness of 10 5 . An externally generated, conical shock wave 
of a strength near that required for separation of the boundary layer was 
used to impose the adverse pressure gradient. Mean- and fluctuating flow data 
are presented for this case. 

The me an- flow data consist of pitot and static pressures and total tem- 
peratures throughout the flow region of interest. In addition, incipient 
separation data obtained from surface flow patterns of alcohol and pressure 
readings from a small orifice dam are presented. These mean-flow data are 
used in conjunction with an existing analytical method to infer the complete 
mean-flow properties including static pressures and flow angles within the 
boundary layer throughout the interaction. The mean-flow measurements serve 
primarily to define the flow in which the fluctuating-flow measurements are 
obtained; however, new information relative to the behavior of the mean flow 
is presented. 

The fluctuating-flow data consist of the turbulent transport terms in the 
conservation equations of mass, momentum, and thermal energy. These quantities 
were determined experimentally with the aid of a hot-wire anemometer using both 
normal and yawed wires. Techniques for the calibration of the wires and the 
interpretation of the hot-wire data are discussed in detail. 

Important conclusions derived from the present investigation are 
summarized below: 

1. The pressure rise to incipient separation is less than previously 
obtained in axially symmetric flow and is much less than corresponding planar, 
two-dimensional flow. 
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2. There are significant static-pressure and flow-angle variations across 
the boundary layer within and downstream of the interaction. 

3. The adverse pressure gradient significantly affects the total 
temperature distribution. 

b. The experimental turbulent stresses agree qualitatively with data 
obtained in incompressible flows with adverse pressure gradients. 

5. The Reynolds normal-stress term normally neglected in boundary-layer 
analyses is important within and just downstream of the interaction. 

6. The turbulent mixing rates are strongly out of equilibrium with the 
mean flow as a result of the interaction. 

7. Because of 6, nonequilibrium modeling of the turbulent mixing rates 
is required, even far downstream of the interaction. 

INTRODUCTION 

The behavior of compressible, turbulent boundary layers in strong adverse 
pressure gradients, particularly those gradients induced by oblique shock 
waves in supersonic flows, is of prime current interest. One major source of 
this interest stems from the desire to predict the behavior of a turbulent 
boundary layer in both internal and external aerodynamic flows on aircraft 
that are to fly at supersonic and hypersonic Mach numbers. Shock waves are 
always present at these Mach numbers; generally, when they interact with the 
boundary layer (whether in an internal flow such as the engine inlet system or 
an external flow such as a deflected control surface), a strong, retarding 

f 

(adverse) pressure gradient acts on the boundary layer. This gradient causes 
a modification of the boundary layer itself, as well as of the flow external 
to the boundary layer. The modifications to the boundary- layer flow through 
the shock-wave - boundary-layer interaction region (figure l) include changes 
in mean profile shape (possibly including separation) and changes in the 
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turbulent mixing properties of the layer. Deviations of the flow external to 
the boundary layer from the classical gas dynamic-discontinuity concept of a 
shock wave include the existence of broad regions of compression and/or 
expansion. For internal aerodynamic flows, the modifications of the external 
flow field are of paramount importance, since the characteristics of the 
reflected compression and expansion waves originating on one wall will deter- 
mine the character of the flow field that interacts with the boundary layer 
on the opposite wall. 

A primary goal of fluid mechanics research today is to develop computing 
techniques that take advantage of recent developments in computer hardware. 

If these computing techniques are to provide realistic solutions for interac- 
tion phenomena, a basic understanding of such phenomena (e.g. , turbulent 
transport rates) is required. The interaction phenomena may be broadly 
divided into two categories denoted here by "mean" and "fluctuating." The 
mean quantities, such as velocity profiles, shock-wave locations, and extent 
of separated regions, are the ones most frequently studied. The fluctuating 
quantities, such as those related to turbulent mixing rates, have been investi- 
gated to a much lesser extent, in spite of their important role in determining 
the mean quantities. 

Existing analytical efforts toward predicting the effects of shock waves 
incident on turbulent boundary layers can be classified into two groups: The 

first considers the interaction phenomenon to be a boundary- layer problem; the 
second considers the interaction phenomenon to be one requiring a separate 
modeling that includes the effects of shock waves. An example of the first 
group is the analysis of Lynes , Nielsen, and Kuhn (ref. l), which is based on 
the turbulent boundary-layer equations. The effects of turbulence are taken 
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into account through an effective eddy viscosity transport model. The agree- 
ment between this method and data obtained in interaction regions becomes 
poorer with increasing streamwise pressure gradients. This may be due to the 
fact that large transverse pressure gradients exist in the interaction region 
(ref. 2), or it may result because the turbulence modeling is not valid for 
these flows. Another example of the first group is the control volume analy- 
sis such as that of Seebaugh et al. (ref. 3) or the improved version by 
Mathews (ref. 4). In this analysis, the boundary layer is assumed to be con- 
fined to a control volume, which is subjected to the conditions of satisfying 
the streamwise momentum balance and continuity relationship. This technique 
has been shown to give reasonable agreement with data for which the pressure 
gradients do not cause large regions of separated flow (ref. 4); however, the 
agreement becomes poorer with increasing pressure gradients. Again, the 
assumption in these methods concerning turbulent mixing and the neglect of 
transverse pressure gradients may cause the poor agreement where pressure 
gradients are strong. Methods in the first group generally yield little 
information concerning the flow external to the boundary layer. 

A highly developed method in the second group is given in reference 5. 
This method assumes that the turbulent boundary layer in the interaction 
region can be divided into two portions: one that responds inviscidly to the 

pressure rise and includes transverse pressure gradients, and one (near the 
wall) that responds as a viscous laminar shear layer with no transverse pres- 
sure gradients. The inviscid portion is computed with a rotational method of 
characteristics program and the viscous portion by a laminar boundary-layer 
program. These two programs are coupled, and the resulting scheme is capable 
of predicting static pressures, total pressures, and flow angles in the 
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interaction region. The method gives results that are in good agreement with 
the observed external flow field and portions of the boundary- layer flow even 
for separated interactions; however, the results derived from the viscous por- 
tion diverge from data when extensive separation is present (ref. 6). Possible 
reasons for this divergence are discussed in reference 6, 

Both the control- volume and the two-layer methods have their advantages 
and disadvantages; they are both ad hoc approximations to the actual physical 
phenomena. These methods have been brought to reasonably "polished" forms 
only within the last three years, and further comparisons with mean-flow data 
are required to assess how good the approximations of each method actually are. 
These two methods treat only modifications to the mean flow and neglect changes 
in the turbulent mixing properties of the boundary layer. However, observa- 
tions (made during the present investigation) of two interactions occurring in 
proximity in the same boundary layer indicate that the modification of the tur- 
bulent mixing properties due to the upstream interaction alters the mean flow 
characteristics of the downstream interaction. This alteration could not be 
accounted for in either of the methods of references k or 5, because they con- 
sider the entire "transfer function" of the interaction to be one of momentum 
removal subject to mass conservation. 

The purpose of the present investigation was to obtain a body of mean- and 
fluctuating-flow data, including the complete Reynolds stress tensor and turbu- 
lent mass- and heat-transfer rates, upstream, within, and downstream of a 
shock-wave - boundary-layer interaction. As a means of determining which of 
the turbulent transport -quantities primarily affect the mean flow, the differ- 
ential conservation equations of mass, momentum, and thermal energy were writ- 
ten in terms of mean and fluctuating quantities and then time averaged. The 
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present work reports the results of an experimental investigation to determine 
each of the turbulent transport quantities and, thereby, to determine which of 
the quantities must be known (or modeled) to predict the changes in the mean 
flow throughout the interaction. 

The measurements were obtained in the turbulent boundary on the wall of 
an axially symmetric nozzle and test section. The nominal free-stream Mach 
number was M = k and the Reynolds number based on boundary-layer thickness 
was Re^ = 10 5 . The adverse pressure gradient was imposed on the boundary 
layer by a conical shock wave generated by a 9° half-angle cone placed in the 
center of the circular test section. The shock strength for this cone was near 
that required for separation of the boundary layer. This experimental arrange- 
ment produces a steep pressure rise in the immediate region of the interaction 
followed by an almost linear, but shallower pressure rise in the downstream 
flow. Mean- and fluctuating- flow data were taken throughout these regions of 
pressure rise. The techniques used for calibration of the wires and interpre- 
tation of both the mean- and fluctuating- flow measurements are discussed. 

REVIEW OF PREVIOUS WORK 
Mean Flow 

Measurements of mean-flow properties of compressible turbulent boundary 
layers in strong adverse pressure gradients induced by shock waves have been 
made by several investigators (refs. 2, U, 7-18). These investigations and 
their primary contributions to the understanding of shock-wave - boundary- layer 
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interaction phenomena have been summarized in references 2, k, lU, and 15 . 

Since mean measurements are made with the intent of improving our understand- 
■'** 

ing of flow in the interaction region, the simplest configurations are chosen 
for investigation. Thus, the investigations of references 2, U, and 7 through 
18 have been performed in two-dimensional flows, either planar or axially sym- 
metric. The planar-flow studies (refs. 2, 7-14) were conducted either on the 
turbulent boundary layer on one wall of a rectangular wind-tunnel test section 
or on the boundary layers developing on flat or curved planar surfaces held in 
a supersonic stream. The studies for axially symmetric flows (refs. 4, 15-18) 
were conducted either on circular wind-tunnel walls or on circular center 
bodies in a supersonic stream. It was tacitly assumed that there was no dif- 
ference in the axially symmetric and planar two-dimensional flows except for a 
coordinate transformation. However, a problem exists in obtaining exactly 
two-dimensional flow on a planar configuration. One such configuration was 
investigated by Lewis (ref. 19) for the case of laminar boundary layers. He 
showed that using planar models with side walls ("boundary-layer fences") did 
not necessarily lead to two-dimensional flow. He also showed that simply 
checking spanwise surface pressure distribution to ensure no large variation 
of pressure across the model at a given longitudinal station is not sufficient 
to guarantee two-dimensional flow. He concluded that an axially symmetric 
model is the only way that an actual two-dimensional flow can be realized. 
Recently, Green (ref. 20) has shown oil- flow photographs of a planar shock- 
wave - turbulent-boundary-layer interaction that indicated a highly three- 
dimensional flow. Similar highly three-dimensional results have been found in 
a supersonic rectangular channel by Reda and Murphy (ref. 2l). The effect of 
these three-dimensional flows on the interpretation of data, presumed to be 
nominally two-dimensional, is unknown. However, it now appears that the only 
practical way to achieve a flow that remains two-dimensional throughout a 
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strong adverse pressure gradient is to use an axially symmetric; configuration. 
Thus, data obtained in axially symmetric environments are given primary con- 
sideration in comparisons with "two-dimensional” theoretical predictions. 

The studies cited above have mapped the behavior of the me an- flow proper- 
ties in shock-wave - boundary-layer interactions as functions of Mach number, 
Reynolds number, and shock-wave strength. However, three areas of interest 
relative to the mean-flow properties remain unsettled: (l) the detection of 

separated flow near a wall; (2) the accurate measurement of static pressures 
in boundary layers with large transverse and streamwise pressure gradients; 
and (3) the accurate measurement of total temperature in such boundary layers. 

The existence of separated flow is usually detected by pitot tube measure- 
ments (as in ref, 8, lU, or 15)* However, a problem with the pitot tube 
arises when it is used near a solid wall. For small separated regions, the 
height of separation can be less than the probe height. For this case, a por- 
tion of the probe is influenced by the impact pressure from the flow external 
to the separation bubble and causes the pitot to read higher than the actual 
pressure that would exist in the absence of the probe. This error causes 
underestimation of the extent of separation. When one is looking for the 
pressure rise to the condition of incipient separation, the rise will be over- 
estimated. To avoid this problem, Reda and Page (ref. 22) used orifice dams 
that in essence determine if the flow next to the wall is moving upstream or 
downstream. 

Static-pressure measurements have been made in shock-wave - boundary- 

layer interactions (e.g. , refs. 2 and 18). However, when static-pressure 

probes are used in regions containing shock waves, the presence of these 

shocks severely affects the probe readings, and as shown in reference 23, the 

magnitude of the error can be nearly as large as the static-pressure change 

across the shock wave. Thus, it would seem highly unlikely that such a probe 
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could accurately detect the pressure gradients in the boundary layer in an 
interaction region. In fact, even the "behavior of the static probe in a 
boundary layer known to have little or no normal pressure gradient (and, of 
course, no shock waves) is not clear. For example, figure 8(a) of reference 2 
indicates a change of approximately 10 percent in pressure from the indicated 
free-stream level as the probe is brought through the boundary layer. In the 
present study, no attempt was made to make direct static-pressure measurements 
in the flow. 

The mean total temperature of a free-stream supersonic flow can be made 
readily by triple-shielded thermocouples. Measurements in boundary layers, 
however, require small probes, and single-shielded or bare-wire probes are 
generally considered. Two types of single-shielded probes have been used. 

The first is the self-aspirated type (ref. 2k or 2) in which the flow through 
the probe, and hence over the thermocouple bead, is governed by the flow con- 
ditions present at the front and rear of the probe. The second is a 
controlled- aspirated type (ref. 25) in which the mass flux past the thermo-- 
couple bead is predetermined by a system of orifices. Since the output of the 
thermocouple depends on its cooling (mass flux past the bead) and on the dif- 
ference between the mean total temperature and the recovery temperature for 
the bead, the advantage of the second type is clear. Both the controlled- 
aspirated probe and some bare-wire probes require a relatively large holder 
and thus may not provide the necessary spatial resolution. 

Fluctuating Flow 

Many investigations *of the fluctuating properties of compressible turbu- 
lent flows have been made (e.g. , refs. lk 9 26-36). The investigations cited 
date back only to about 1950 when electrical compensation networks were 
applied to constant-current anemometers to make them a serious research tool 
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in high-speed turbulence. Kovasznay (ref. 37) summarizes the development of 
constant-current equipment leading to frequency responses in excess of 50 kHz. 
Constant-temperature anemometers that at least double or triple this frequency 
range are now more or less "off-the-shelf" items. Methods other than the hot- 
wire anemometer for making fluctuating measurements - for example, the 
electron-beam technique studied by Wallace (ref. 32) - have been investigated 
recently. However, Kovasznay (ref. 38) concludes that the hot-wire anemometer 
is, at present, the best instrument for this purpose, and it was used in the 
present study. 

The now classic works of Kovasznay (refs. 28 and 30 ) and Morkovin (ref. 
2J) led to the relatively sophisticated present status of signal interpreta- 
tion in high-speed flow (summarized by Morkovin, in ref. 39)* All the cited 
references use these hot-wire signal-interpretation techniques to separate the 
fluctuations in mass flux and total temperature. The wire is operated at at 
least three different wire temperatures (so-called "over-heat" ratios). The 
sensitivity of the wire to the various fluctuations changes with the over-heat 
ratio, and the wire thus may be used to obtain the desired fluctuations in an 
unknown field. In addition to the normally predominant entropy and vorticity 
fluctuations in a boundary layer, there are pressure fluctuations ("sound" 
mode or "aerodynamic noise") that usually become significant external to a 
viscous boundary layer. Investigations of this sound mode and its interaction 
with shock waves and the other modes have been investigated (refs. lb 9 27, 33, 
36, UO-UU). Again, the techniques of reference 27, 28, or 30 are used to 
separate the sound from the entropy and vorticity. These techniques have been 
in existence for over 15 years and are not discussed here in detail; Grande 
(ref. lU) gives an excellent review of this work. 
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The results of these investigations are reviewed here as they apply to 
the determination of fluctuations in the turbulent boundary layer and near- 
field external flow. 

Kistler (ref. 29) has shown distributions of <(pu)’> } <T!>, and (pu) 1 T? 

"C b 

across a zero-pressure-gradient turbulent boundary layer at Mach numbers of 
1.71, 3.56, and U .67 for a single streamwise location. Laufer (ref. 36) has 
measured the mass-flux and total-temperature fluctuations external to a zero- 
pressure-gradient turbulent boundary layer and concluded that they were pro- 
duced solely by pressure fluctuations radiated from it. These two investiga- 
tions were carried out using a single, normal wire. However, none of the 
measurements was converted to a form that could be used to assess relative 
magnitudes of the turbulent stresses, such as pu’ 2 or up f u T , in the momentum 
equations. Morkovin and Phinney (ref. 26) considered the use of a yawed wire 
to determine the v’ 2 -like fluctuations and the shear-stress-like term u’v'. 
They present the shear-stress measurement at one point in a zero-pressure- 
gradient turbulent boundary layer. After ik years, this apparently is still 
the only data point of its kind. 

Grande (ref. lb) used a single normal wire to obtain fluctuating flow 
measurements in the turbulent boundary layer and external flow in an oblique 
shock-wave (turning angles of 6°, 8°, 10°, and 12°) - turbulent -boundary-layer 
interaction at a Mach number of 2.5. This is certainly the most detailed 
survey of the most general flow configuration studied to date. He presented 
values of <(pu) f >, <t>, <a>, and <tt>, and the fluctuating pitot pressure 

obtained with a Kulite sensor, at several points throughout the inter- 
action. However, most of the data were obtained outside the boundary layer. 

He showed that the fluctuations downstream of an interaction were signifi- 
cantly changed due to the interaction. These findings are summarized 

schematically in figure 100 of reference lb. 
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Grande used an anemometer system whose upper frequency response was at 
most bo kHz and, therefore, his measurements only encompass fluctuations whose 
scales are larger than about 6. This system is probably inadequate for 
studies of boundary-layer turbulence but may have been adequate for his pur- 
poses. In spite of the limitations of Grande’s system, he was able to show 
that existing signal-interpretation techniques for hot-wire anemometers using 
normal wires could be used with assurance in very complicated flow fields. 

As noted previously, the turbulent mixing properties of compressible, 
turbulent boundary layers have been investigated to a much lesser extent than 
the mean- flow properties. Furthermore, turbulent shear stresses, which could 
likely play an important role in altering the boundary- layer flow downstream 
of an interaction, have never been experimentally investigated. 

ANALYTICAL CONSIDERATIONS 

We will consider the flow of air (a compressible, viscous, heat- 
conducting medium) at pressures and temperatures for which the continuum flow 
assumptions are valid and no dissociation of the gas takes place. In addition, 
we assume the gas behaves as a Newtonian fluid (i.e., a linear stress - rate- 
of-strain relationship is an adequate assumption for the flow process involved) 
and that Stokes’ hypothesis is valid for these processes (y \ = ~2/3u; i.e., 
the relaxation times for energy transfer from translational to rotational or 
vibrational modes in the molecules are essentially zero). Further, we con- 
sider only flow without body forces. We will consider the equations of motion 
in cylindrical coordinates and make the assumption of axially symmetric flow. 

Of course, these equations reduce to the equations of planar, two-dimensional 
flow in the limit as r -* 00 . The coordinate-system nomenclature used through- 
out this work is shown in figure 2. The axial coordinate is denoted by x; 

the radial, r; and circumferential, 0. The velocities in these directions are 
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denoted "by u, v, and w, respectively. In addition, y denotes the direction 
negative to that of r and is measured from the wind-tunnel wall rather than 
the centerline. 

Using these coordinates, the governing equations, which apply instanta- 
neously to the local flow field, may he written for the flow as follows: 
Continuity : 

ifi. + + 1 9 . r P. v . = o 

at 3 x r 3 r 

Momentum : 
x: 

3 3 i* x x 

3 pv Bpvu l_ 8 rpw pww _ xr + 1 _ rr BO 
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We now examine explicitly the "turbulence" and use the Reynolds 
substitutions 

u=u+u* p=p+p’ h=h+h* 
v = v + v ! p = p + p 1 y = y + p ! 
v = w + w ! T = T + T 1 $=$+$» 


where the time average of any fluctuating quantity is zero (e.g., u f = 0). 

When these quantities are substituted in the previous equations, the time 
average taken, and the flow under consideration is assumed to be steady in the 
mean (i.e., 3pu/3t, etc.= 0) with w = 0 for axially symmetric, two-dimensional 
flow, we obtain the following system: 

Continuity : 


Bpu .L 3rpv + 3p ? u ! + 1 Brp* v' 
Bx r Br Bx r Br 


(i) 


Momentum: 


x: 


Bpuu + 1_ Brpuv _ _Bp_ + _3r 

Bx r Br Bx Bx T xx 


(pu’u f + 2up f u* ) ] 


+ — —■ r[x - (pu’v 1 + up ? v T + vp'u r )] (2a) 

r Br rx 


Bpuv 1 3rpw _ 
Bx r 3r 


— [ T — oW * W 1 1 

Br r LT ee pW J 


-- [t -(pu 1 v f +up T v ? +vp ? u* )]+ — ~ r[x -(pv f v’+2vp 1 v* )] (2b) 

Bx xr r Br rr 


15 


6: 0 = 9x ^ T xe “ ^ pu ' w ' + up'w')] + ^7 97 r2 [ T r e “ (p v,w ' + vp'w' )] 


Energy : 


)£uh + 1 9r gvh = - 3 £ + -|E +u , i£l +v , i£l + ^ 
3x r 3r 3x 3r 3x 3r M 


*3x ^ fx " ( pu?hf + u P ?hf + h P* 

1 _3 ~ 3T 

r 3r r 3r 


u») 


(pv’h T + vp ! h f + hp'v’ 




State: 


(2c) 


(3) 


p = H(pT + p ? T 1 ) 


w 


where the t. , k( 3T/3x) , and k( 3T/3r) terms still contain correlations 
between y 1 and 3u ! /3x-like terms and the k ! and 3T ! /9x-like terms. 

We now consider some of the individual terms in equations (2) through (4) 
together with experimentally observed features of supersonic flows to reduce 
(slightly) the complexity of this governing system of equations. First, con- 
sider the following definitions: 

T = t - (pu'u 1 + 2up T u ? ) 
xx xx 

T = t - (pv f v f + 2vp f v’ ) 
rr rr 

T = T - (pu f v* + up'v 1 + vp f u* ) 
rx rx 


T 


00 



pw T w f 


Now, if the stress tensor is to remain symmetric, we must have T = T . We 

3 t/ ^ xr rx 

also have 


T ex = T X0 = T xe - (pu ' w ’ + ' Jp ' w ' ) 


T A = T ^ = T A - (pv f w f + Vp f w ? ) 
0r r0 r0 


Thus (2c) may be written: 


o = + — -1 r 2 T 

u 3x r 2 3r r r0 
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Now, for an axially symmetric flow (e.g. , ref. 1+5), or, for that matter, a 
planar, two-dimensional flow (e.g., ref. 26), T must he identically zero. 
We then have: 0 = (3/3r)r 2 T r0 . Integrating from r = 0 to r we see: 


0 - f 


r=o 


a? r 2 T re dr = r2T re - r2T re 


+ f(x) 


r=o 


Assuming the physical variables in T A to he hounded, we have r 2 T I = 0 

r0 r0 I r=o 

and thus must also he identically zero. With this result, the 9 momen- 

tum equation is degenerate. 

Thus , we may write the equations describing axially symmetric, turbulent 
flow in the following form: 

Continuity : equation (l) 

Momentum: 


x: 


r : 


6: 


9puu 1, Brpuv _ _3 jd 

9x r 3r 3x + 



1 

r 


3rT 


xr 

3r 


3puv + 1 3rpvv _ 3j3 3T rx 1 3rT rr T 00 

3x r 3r 3r 3x + r 3r " r 


r0 


- T X0 = 0 


(5a) 

(5b) 

(5c) 


Energy : equation (3) 

State : equation (1+) 


These equations are exact in that no order-of-magnitude approximations or any 
other assumptions relative to the size of the terms have been invoked. We now 
make use of existing knowledge of the size of some of these terms to reduce 
the complexity of equations (3) and (5), and thereby see which terms in these 
equations are significant in altering the mean flow. First consider the term 


where 


T = U 
XX 


2 £-!<*•*) 


+ p' 


2 ^-§<v • ?)• 


V • v 


3u ( 1 3rv 
3x r 3r 


IT 


and 


(V 




1 3rv f 
r 3r 


To simply neglect the second part of x is inconsistent with a nonzero h f , 
since y f is a (weak) function of h f . However, indications are (see Laufer, 
ref. 46, or Favre, ref. 47) that gradients of the fluctuating viscosity- 
velocity gradient correlations are usually much smaller than the gradients in 
the mean terms so that in all the x terms we need only retain the terms 
involving the mean viscosity and velocities. Thus, we write the explicit 
relations : 


T = y 
rx 


X0 


r0 


_ 2 3*' 

■ 3 

f? 

2 

L 2 ^ 

“ 3 

r 2 i 

2 

L 2 ^- 

3 

ii + 

3v 

[aX 

3x_ 

(pu ! w f 

+ 

(pv’w 1 

" + 


- (pu’u 1 + 2up f u ? ) 


(V • V) 


- (pv T v ! + 2vp ’ v 1 ) 


V) 


- pw’w 1 


- (pu’v* + up’v 1 + vp’u 1 ) 


(6a) 
(6b) 
(6c) 
(6d) 
(6e ) 
(6f) 


Now, considering the energy equation (3), we first examine the nature of 
the terms 


u 



and 



Based on measurements by Kistler (ref. 29) and the discussion by Laufer (ref. 
46) for flows below a Mach number of 5, one would expect that u ! ( 3p ? /3xT and 
v f (9p f /3r) are small compared to their mean counterparts u(3p/3x) and 
v(3p/3r), and may be neglected. Next we write the terms: 


k 



31 + k , ill 

3x 3x 
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and 



3T 

3r 


+ k f 


3T 

3r 


Again, based on arguments similar to those that lead to neglecting U' (9u'/9x), 
etc., we assume that k'(9T'/9x) and k*(3T'/9r) are small. Finally we con- 
sider the term y$. Recalling the definition of y4> as the tensor product 


we may write 


= t : VV 


8u. 

TiJ ^ 





3u. 

1 



+ 


ij 


9u' 

1 



contain fluctuating terms such as y' and 3u'/9x, which are 
in the term r| ” ( 9u!y 3x^ ) . This term contains double and 
triple correlations of the viscosity and velocity gradients , which we have 
previously neglected. Therefore, it is not inconsistent to neglect the term 
T iJ^ 9u i/ 9x ,j^ and sim P 1 y write y? = y$. In view of the foregoing argument, 
equation (3) may be cast in the following form: 


The terms t! . 

ij 

then correlated 


where 


and 


3puh + 1 3rpvh 
9x r 9r 




3p - 

+ v 


3x 


(pu’h 1 


(pv ! h T 




+ up , h f + hp f u f ) 
+ vp f h* + hp^ v" 1 ) 


9rQ r 

3r 


( 7 ) 


We now consider the terms such as pu f v ! and pv'h' that contain triple 
correlations; for example, 

pu f v f = pU*V f + p’u’v* 

Laufer (ref. h6) gives reasonable arguments that lead to 
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< ~5 percent 

pu 1 V* 

In addition. Way and Libby (ref. 31 ) have recently reported measurements of 
the triple correlation p'u'u* in a low- velocity flow that contained signifi- 
cant density fluctuations. Their reported data indicate 

p T u f u T < 5 percent 
pu f u f 

Considering the expected accuracy to which quantities such as pu 1 ^ or 
pu f v f ~ can be measured, it thus seems reasonable to neglect the triple 
correlations. We had previously considered the correlations 




u 


8x 


and 



to be small. For weak (either unseparated or weakly separated) shock-wave - 
boundary-layer interactions, such as considered in this study, Grande (ref. 
llO has shown that the pressure fluctuations themselves induced by the inter- 
action are small in comparison with the boundary- layer turbulence. Now we 
examine the consequence of assuming the fluctuations p 1 to be zero. This 
assumption is also an important basis for the hot-wire anemometer data reduc- 
tion procedure discussed later. For p f = 0, it can be shown by logarith- 
mically differentiating the equation of state: 



hi 

h 


this is used in the following. The effect of neglecting the triple correla- 
tions and pressure fluctuations on the turbulent heat-transfer terms is dis- 
cussed below. First, neglecting triple correlations, we have: 


Q = k - ( pu T h* + up 1 h f + hp'u 1 ) 


(8a) 
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(8b) 


% = k - (pv'h f + vp f h f + hp 1 v T ) 


Consider, for example, the first term inside the parentheses of and employ 

the relationship p’/p = -h f /b based on the assumption of zero pressure 
fluctuation. We have: 


and similarly: 


pu»h T = -pu'p f =r 

P 


-hp f u ! 


pv T h f = -hp'v* 


Thus we see that, for the form of the convective terms (left-hand side) used 
in equation (T), the expressions for the heat-transfer rates reduce to: 



up’h' 


(9a) 


Q 


r 



vp’h' 


(9b) 


Making use of the above assumptions, we obtain the conservation equations 
and the equation of state in the forms given below: 


Continuity : 

3pu . 1 9rpv a~p v u y , 1 ar^'v' _ n 
r 3r 9x r 3r 

Momentum : 

x: 9 PU 2 * 1 9rpuv _ 3^ 3T xx 1 3rT rx 

r 3r 3x 3x r 3r 

r: 9puv ( 1 3rpv 2 3p 3T rx 1 3rT rr T 66 

3x r 3r 3r 3x r 3r r 


( 10 ) 


(lla) 

(lib) 


6: 

where 


T = T = 0 
r0 x0 


T = y 

XX 


T = u 
rr 


T ee ^ 


2 £-!<* •*> 


2g-f (V . *) 


- (pu f 2 + 2up f u f ) 

- (pv 1 2 + 2vp f v f ) 


2 5 - f (7 • 5) 


- (pw' 2 ) 


(11c) 


(12a) 

(12b) 

(12c) 
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(I2d) 


rx 


- 3v 

= W [^7 + 


- (pu’v* + up’v 1 + Vp’u') 


T „ = -(pu'w ! + up f w ! ) 
x0 


Energy : 


= -(pv’w f + vp'w 1 ) 


— — 3Q 3rQ 

+ 1 3r fi. v k = u^+v-^+y $ + -rr L + - r 

3x r 3r U 3x V 3r y * ^ - 


8x r 3r 


where Q and Q are given by equation (9). 
x r 


State: 


= R(pT + p'T' ) = RpTll - 


(I2e) 
( 12f ) 

(13) 

(lU) 


An experimental investigation was conducted to measure the mean flow 
properties and each of the turbulent transport terms in equations (10) through 
(lU) in a shock-wave - boundary-layer interaction. The remainder of the thesis 
reports the results of that investigation. 


EXPERIMENTAL INVESTIGATION 


The experimental investigation was conducted in an axially symmetric flow 
facility to minimize three-dimensional effects. Figure 3 is a schematic dia- 
gram of the experimental facility used in the investigation. A complete des- 
cription of the experimental apparatus is given in appendix A. 

The turbulent boundary layer under investigation was the one that devel- 
oped on the nozzle and test section wall. The boundary layer upstream of the 
interaction is referred to as the "undisturbed" boundary layer. At the begin- 
ning of the interaction it had a history of both a favorable and a slight 
adverse pressure gradient (appendix A). The Mach number for the core flow in 
the test section was 3.88 ± 0.02. The Reynolds number based on the thickness 
of the undisturbed boundary layer was 8.7 x lO 1 *. The total temperature was 
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5l+0° R (300° K), and the wall temperature was near that of the adiabatic wall 
temperature for the test section Mach number. 

The boundary layer was subjected to an adverse pressure gradient generated 
by a cone placed on the centerline of the test section. The resulting wall 
static pressure distribution is depicted in figure 3. Three rather distinct 
pressure gradients are evident: the nominally zero pressure gradient ahead of 

the incident shock wave, and two adverse pressure gradients, both of which are 
studied in the present investigation. The second gradient is a rather steep 
one induced directly by the incident-reflected shock-wave system, while the 
third (downstream of the second) is a less severe, nearly constant gradient 
induced by the conical flow field. 

The magnitude of the pressure gradients can be varied by using different 
half-angle, straight cones. It was decided that a flow for which the imposed 
pressure gradient was near that for incipient separation should be investi- 
gated. As a preliminary step in the experimental investigation, then, the 
pressure gradient required to just separate the boundary-layer flow had to be 
determined. This was done by a technique using surface flow patterns of alco- 
hol in conjunction with pressure measurements from a miniature orifice dam 
(appendix B). With the use of this more sensitive technique, the incipient 
separation pressure rise was found to be much lower than that reported for com- 
parable Mach and Reynolds numbers in reference 15. The results of the incipient 
separation study showed that a cone of 9° half-angle produced a pressure gradi- 
ent that was very near that required for incipient separation, so the 9° cone 
was used in the present study. 

The results of the experimental investigation consist of mean- and 
fluctuating-flow data obtained throughout the boundary layer and external flow 
at the locations shown in figure The mean-flow data consist of pitot and 
surface-static pressures and total temperatures. The pitot pressures were 
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obtained with a pitot probe, and the total temperatures were obtained by using 
a hot wire as a resistance thermometer. The fluctuating- flow data consist of 
the turbulent transport terms for mass, momentum, and thermal energy. The tur- 
bulent terms were obtained by using a constant-temperature hot-wire anemometer. 
Both normal and yawed wires were used. A detailed discussion of the instru- 
mentation is given in appendix A. 

The experimental operating procedures and data— reduction techniques for 
both mean- and fluctuating- flow measurements are given in appendix C. Special 
attention is given to describing the theory and interpretation of the hot-wire- 
anemometer data in appendix D, and the data accuracy is discussed in appendix 

E. 


RESULTS AND DISCUSSION 

The principal experimental results are shown in figures 5 through 17 • 

The data are also presented in tabular form to facilitate the use by others 
who may wish to examine aspects of the data which are different from those dis- 
cussed below. The mean-flow data are listed in table 1 and the fluctuating- 
flow data in table 2. 


Mean-Flow Data 

Mean- flow data were obtained at a single Mach number and Reynolds number 
(M = 3.88; Re 6 = 8.7 x 10 4 ) and a single shock-wave strength (cone half-angle 
of 9°). The pitot pressures and total temperatures obtained through the 
region of interest (fig. 4) are shown in figure 5. The incident shock wave is 
evident in plots of the pitot pressures given in figure 5(b) through (d), and 
the reflected shock-wave system is evident in figure 5(g) through (j). The 
temperature profiles have the shape expected for an adiabatic flow. The tem- 
peratures are lower near the wall than in the free stream and higher than the 
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free-stream value in the outer part of the boundary layer. These profiles are 
qualitatively the same as those found in reference U8. The wall temperature 
is the adiabatic wall temperature for M = 3.88, assuming a recovery factor of 
0.90. The height at which the temperature becomes essentially that of the free 
stream is denoted by 6 T-t , the boundary- layer thickness based on the total- 
temperature profile. The pitot-pressure profiles reveal a well-defined 
boundary— layer tnickness shead of the incident shock wave. Downstream 

of the incident wave, however, the Mach number in the inviscid flow is no 
longer uniform, so that a boundary— layer thickness based on Mach number pro- 
files is not reliable. Thus, the term "boundary— layer thickness" used in the 
present study refers to 6 Tt , which agrees with the 6 pt ^ upstream of the 
incident shock and which remains clearly definable throughout the interaction 
and downstream flow. 

The data presented in figure 5 can be combined to give the representation 
of the mean-flow field shown in figure 6. The incident shock wave is shown as 
a single line (discontinuity in the flow) even though from pitot— pressure pro- 
files alone the shock appears as a band. The hot-wire anemometer was used, as 
discussed in appendix C, to locate the shock more accurately. The line repre- 
senting the incident shock wave has been drawn through the experimental points 
obtained by the hot-wire technique. These points are located approximately in 
the center of the band given by the pitot-pressure survey. The reflected 
shock-wave system is shown as two separate compression bands. The presence of 
two compression bands could not be shown conclusively based on the pitot- 
pressure profiles alone; however, with the use of the hot-wire anemometer it 
was found that two bands were present. Since the shock wave used to generate 
the pressure rise is believed to have been insufficiently strong to separate 
the boundary layer, the term "separation shock wave" probably should not be 
used to describe the upstream (first) reflected compression band. (Even if 
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there were a small separation present, the we 11- developed, two-compression 
configuration could probably not be explained on the simple basis of the exis- 
tence of a region of separation. ) This upstream band is denoted by the term 
"induced shock wave," indicating that it is probably induced from the turning 
of the supersonic flow within the boundary layer due to a thickening of the 
subsonic portion of the layer caused by the forward feeding of the pressure 
rise. Also indicated in the flow field is the edge of the boundary layer 
determined from total-temperature measurements. The induced wave merges with 
the reflected wave soon after leaving the boundary layer. The boundary-layer 
thickness is reduced by about 35 percent in passing through the incident- 
reflected shock system. The thickness of the layer downstream of the reflected 
shock remains constant even though the pressure is rising. 

The experimental surface-static-pressure distribution is shown in figure 
7. A region of nominally zero-pressure-gradient flow exists upstream of the 
incident shock. This region is followed by a steep pressure rise caused by 
the incident-reflected shock system and then by a shallower rise associated 
with the conical flow field. No attempt was made to measure static pressures 
in the flow field because of the large uncertainty present in interpreting 
results of static-pressure-probe data when shock waves are present (see Review 
of Previous Work). However, a technique is outlined below that allowed 
reasonable estimates of the static pressure throughout the flow field. 

The static pressure can be computed from the Rayleigh pitot formula using 
the measured pitot pressure if the actual total pressure is known. Within the 
boundary layer, the total pressure is what we seek to find and is, of course, 
unknown. However, outside the boundary layer the total pressure may be known 
(e.g. , from inviscid-flow calculations). In the particular case of a shock 
wave generated by a 9° half-angle cone, the total pressure behind the incident 
shock is 99.6 percent of the tunnel stagnation pressure. The reflected 
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compression system consists of two broad compression bands (and not a single 
shock wave) so that there is probably very little further loss in total pres- 
sure. Therefore, it is reasonable to estimate that the total pressure loss 
across the shock system is at most 0.5 percent. (Even if the incident shock 
were reflected as a single shock, there would be a net loss in total pressure 
of only 0.7 percent.) In view of the almost negligible losses across the 
shock systems, the static pressures in the flow field external to the boundary 
layer were computed from the Rayleigh pitot formula assuming no total pressure 
loss. These pressures are shown in figure 8, where they have been normalized 
with respect to the measured value of the local wall static pressure. Also 
shown in figure 8 are the static pressures as computed by the method of refer- 
ence 5, which was discussed in the Introduction, The calculation procedure of 
5 was started at station 2.60 in. ( 6.60 cm) (fig. 8(b)) where 
p / p wall was taken "k 0 be 1*00 across the boundary layer. The static pressures 
determined in this manner agree to within 2 percent of those computed from 
pitot pressures except in the vicinity of the reflected compression system 
(e.g., fig. 8(g)). Here, the method of reference 5 predicts a shock wave 
rather than the broad compression indicated by the experimental data. It can 
be seen from figure 8 that there are large static pressure gradients normal to 
the flow within the boundary layer. In figure 8(z) , the expansion from the aft 
corner of the shock-wave generator is evident. At this station (x = k.60 in. 
(11.68 cm)), the expansion is still outside the boundary layer so that the data 
discussed should be free of expansion effects. 

The predictions for static pressure below 6^ seem to be reasonable; 
however, a check on the credibility of these predictions is needed. To this 
end, the method of reference 5 was used to predict pitot pressures. These were 
then compared directly with the pitot pressures measured in the interaction 
region below 6^ (fig. 9). The starting profile at x = 2.60 in. (6.60 cm) 


27 



is shown in figure 9(a). Experimental data were used for points within the 
boundary layer and below the shock wave, while points above the wave were 
obtained from a theoretical conical flow solution. Profiles are also shown for 
X = 3.00 in. (7.62 cm), 3.20 in. (8.13 cm), and 3.60 in. (9.14 cm) (figs. 9 (*b ) 
through (d)). At x = 3.60 in. (9.lU cm), the reflected compression is just 
outside the boundary layer. The predictions agree quite well with the data, 
with the largest difference of about 10 percent observed to occur near 
at x = 3.60 in. (9.lU cm) (fig. 9(d)); in most instances, the predictions are 
within 5 percent of the measured values. In view of the agreement between 
predicted and measured values of pressure, it was decided that the method of 
reference 5 provided reasonable estimates of the static pressure below 
These values have been used for points within the boundary layer, while at the 
wall and outside the boundary layer the experimentally determined values have 
been used. 

Estimates of the flow angles throughout the flow field were also obtained 
from the method of reference 5. Flow angles must be input at the initial sta- 
tion; however, no attempt was made to obtain them by measurement. Above the 
incident shock wave, the flow angles were taken to be those determined from the 
theoretical conical flow solution. At (Srp^ the flow angle was estimated from 
the continuity equation by assuming planar flow in a zero pressure gradient and 
integrating the continuity equation to give v e / u e = ^<5*/dx. The flow angle 
was then assumed to vary linearly from the value at 6^ to zero at the wall. 
Experimental values of 6* in the undisturbed boundary- layer flow are shown 
in figure 10. The rate of change of 6* indicates a flow angle of about -0.4 
at (Sip^. Using this initial set of flow angles, the method of reference 5 
predicts the angles throughout the flow field shown in figure 11. As an 
approximate check of the predicted values, estimates were made of flow direc- 
tion at the boundary-layer edge in downstream flow (fig. ll(g) through (&)) "by 
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use of the continuity equation. Since, as noted previously, the houndary- 
layer thickness does not change downstream of station 3.60 in. (9.1*+ cm) (fig. 
11(g)), the flow angle at can be estimated by using the values of the 

mass flow entrained into the boundary layer. The mass flow in the boundary 
layer, normalized with respect to the undisturbed boundary- layer flow 
(x = 2. 40 in. (6.10 cm)), is shown in figure 12. The flow angles determined 
in this manner are shown in figure ll(g) through (&). The agreement between 
the predicted flow angle and the angle determined from mass entrainment is 
quite good at x = 3.60 in. (9.1*+ cm). Further downstream, however, the pre- 
dicted values become successively larger than those based on entrainment. 

Since the angles determined by entrainment are only estimates (probably accu- 
rate to within only 1°, if that) it was decided that the predicted values 
should be used for the mean radial velocity v. As can be seen in figure 11, 
flow angles within the boundary layer may be quite large. 

The total-temperature data can also be used to examine the effect of a 
strong adverse pressure gradient on the temperature profile. To this end, the 
temperature data are plotted in figure 13 on the nondimensional temperature 
versus velocity coordinate. The linear relationship for zero-pressure-gradient, 
unity-Prandtl-number flow (the so-called Crocco relationship) is shown along 
with the data for four streamwise stations. The entering profile (x = 2.60 in. 
( 6.60 cm)) has a behavior typical of nozzle-wall (favorable pressure gradient) 
temperature distributions (see ref. *+8). However, a rather pronounced change 
in the profile occurs after the flow encounters a region of adverse pressure 
gradient (x = 3.20 in. (8.13 cm)). The profiles continue to change as the 
flow moves into regions of higher and higher pressure. Thus, the departure of 
the temperature from the Crocco relationship is substantial, just as it is for 
flows undergoing sustained favorable gradients (e.g., the nozzle wall). 
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One final topic relative to the mean flow can be discussed with the aid of 
figure 12. The analytical method described in reference U, for example, assumes 
that the mass addition to the boundary layer throughout the incident-reflected 
shock system is negligible. The data of the present study indicate that this 
condition is true as can be seen from figure 12. The rate of mass addition up 
to station x = 3.60 in. (9.1*+ cm) is essentially the same as in the undisturbed 
boundary layer; however, downstream of this station the rate of mass addition is 
more than an order of magnitude larger than that in the undisturbed flow. 

Fluctuating-Flow Data 

The fluctuating- flow data consist of the turbulent mass-transfer , 
momentum-transfer (Reynolds stress), and heat-transfer terms discussed in rela- 
tion to equations (8) and (10) through (13). Only the data relative to the 
continuity and momentum equations are discussed here. However, all the 
fluctuating- flow data (including the heat-transfer term pv ? T* , which appears 
in eq. (8), but not eq. (9)) are presented in table 2. We examine the continu- 
ity and momentum equations to determine which of the fluctuating terms are sig- 
nificant in altering the mean flow throughout the interaction region. 

The turbulent mass-transfer terms p’u’ and p T v f in the continuity equa- 
tion (10) are shown in figure lU. The values of p f v* actually go to zero at 
6 ij , while those of p'u* have a small but positive value in the free stream. 
The wall values are assumed to be zero, consistent with p’, u* , and v* each 
being zero at the wall. The shape of the p*v f profiles remains qualitatively 
the same throughout the interaction and in the downstream flow, attaining 
maxima near 6 T ^/2. The profile of p T u', however, changes qualitatively 
throughout the interaction. The maximum exists near the wall in the upstream 
(zero pressure gradient) flow. Within the interaction, this maximum near the 
wall is washed out. However, the downstream flow appears to be relaxing to a 
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profile (with a maximum near the wall) similar to that of the upstream flow* 

In addition to the observable qualitative changes, significant quantita- 
tive changes are evident in figure lb m There are significant increases in 
both p’u 1 and p'v 1 downstream of the incident shock and within the boundary 
layer (fig. lb (d)). Further increases are evident downstream of the reflected 
shock system. These increases persist in the downstream flow to the last mea- 
suring station (x = H.60 in. (11.68 cm)). The slight increases in the magni- 
tudes of p f u ! and p'v 1 evident near the reflected shock system (fig, lb(g) 
and (h)) are probably due to pressure fluctuations associated with the 
reflected shock system. These pressure fluctuations are neglected in the 
present study as discussed in the Analytical Considerations Section. We now 
turn directly to the derivatives of p f u* and p'v 1 in equation (lO). Figure 
15 shows the experimental values of Sp’u’/^x and (l/r ) ( 3rp f v* /3r) at four 
streamwise stations. Figure 15(a) shows data at an upstream station (x = 2.b0 
in. (6.10 cm)), 15(h) at a station within the interaction (x = 3.20 in. (8.13 
cm)), 15(c) at a station just downstream (x = 3.80 in. (9.65 cm)), and 15(d) 
at a station far downstream (x = U.60 in. (11.68 cm)). The derivatives were 
obtained from the data shown in figure lb by first smoothing the data and then 
differentiating the smoothed data as discussed in appendix C. 

It is clear from figure 15 that the term 3p T u f /3x is almost everywhere 
at least one order of magnitude less than (l/r) ( 3rp T v T /8r ) . Furthermore, the 
difference between the two derivatives becomes more pronounced as the flow 
goes through the adverse pressure gradient. Thus, for flows similar to the 
one in the present study, it appears that the term 3p ? u f /3x may safely be 
neglected in analytical descriptions of such flows. Invoking this assumption, 
equation (10) may be written as: 

3pu + 1 3rpv + 1 3rp 1 v 1 = 

3x r 9r r 3r 
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or _ 

3£u 1 9££V = 

Sx r 

which is just the usual continuity equation obtained hy employing the boundary- 
layer approximation. 

All the turbulent moment urn-transfer terms (Reynolds stresses) that appear 

in equation (ll) except were measured in the present study and are shown 

in figure l6. All values of T , except T at the wall, are taken to be 

x j rx 

just the turbulent components of the stresses. This was justified by calculat- 
ing the molecular contribution to from the mean-flow velocity data. At 

measuring stations away from the wall, the molecular portion of the total shear 
stress was found to be less than the turbulent portion by at least two orders 

of magnitude. At the wall, all the T. . are assumed to be zero, except T . 

i j x x 

The wall values for T were calculated from the value of c derived from 

rx f 

a fit of the mean- flow velocity profile to a compressible form of the law of 
the wall (ref. k ) . 

The qualitative behavior of these profiles can be discussed in light of 

known results obtained for incompressible, adverse-pressure- gradient flows 

(e.g. , Sandborn and Slogar, ref. 50, or Bradshaw, ref. 51). 

The general behavior of the T profiles does not significantly change 

rx 

throughout the interaction. The values go to zero at the edge of the boundary 
layer, have a maximum near 6/2, and have approximately constant values from 
the wall to the first data point away from the wall. However, the magnitude 
of T x increases markedly as the flow goes through the interaction. This is 
in qualitative agreement with results shown in references 50 and 51 for 
incompressible flows in adverse pressure gradients. The behavior of T^ x in 
the undisturbed flow (e.g., x = 2.40 in. (6.10 cm), fig. 16 (a)) is similar to 
that shown in references 50 and 51 for a mild adverse pressure gradient. As 


32 


noted in appendix A, the present "undisturbed" flow has, in fact, gone 
through a slight adverse gradient. 

The behavior of T^, T^, and T 00 are also similar to their incompres- 
sible analogs. The values of T are larger than either T or T in the 

xx rr 00 

lower part of the layer; T is less than T near the wall, while T is 

greater than T 00 in the outer portion of the boundary layer. 

The value of T 0x was measured in the present study. As was shown for a 
flow that is exactly two-dimensional, this term must be zero. It can be seen 
that the values of T 0 ^ are at least one order of magnitude lower than those 

of T rx over nearly the entire flow. Although the values of T 0x are small, 

they are not zero. Furthermore, the values of T 0x exhibit a systematic pat- 
tern at all the measuring stations. The values are positive in the portion of 
the layer near the wall, and negative in the outer portion. The pattern is 
accentuated as the flow passes through the adverse pressure gradient (and as 
the boundary layer becomes thinner). This behavior might possibly be explained 
by the presence of streamwise (Goertler) vortices that are approximately the 
scale of the boundary-layer thickness. No definite conclusion relative to this 
possibility could be drawn from the present data and further work is required 
before the cause of these patterns can be confirmed. One can conclude, how- 
ever, that the flow remained nominally two-dimensional throughout the adverse 
pressure gradient. 

With a view toward determining which of the stresses play a significant 
role in altering the mean flow in equation (ll), we turn directly to the gradi- 
ents of the stresses. The gradients of mean pressure also appear in the 
momentum equations. Thus, to determine which gradients alter the mean flow, it 
is necessary to compare the turbulent stress gradients with themselves and with 
the pressure gradients. 
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The derivatives that appear on the right-hand side of the x-momentum 
equation (eq. (lla)) are shown in figure IT. These derivatives were obtained 
from the data shown in figure 1 6 by the same technique used to obtain the 
derivatives shown in figure 15* The same four stations as shown in figure 15 
are used to show representative data for the entire flow* These are upstream 
(fig. 17(a)), within (fig. 17(b)), just downstream (fig. 17(c)), and far down- 
stream (fig. 17(d)) of the interaction. 

In the upstream, or undisturbed, flow (fig. 17(a)) 3p/3x is nominally 
zero, while 3T /3x is about 10 percent of (l/r)(3rT /3r). Thus, the 

9 XX rJt 

upstream flow conforms reasonably well to the boundary- layer approximation (in 
which 9T xx / 9x is assiime{i to be zero)* The behavior is sharply altered 
within the interaction region as is evident in figure 17(b), where the value 
of 3p/3x is larger in magnitude than either turbulent stress term in the 
lower half of the boundary layer, while in the upper half dp/dx and 
(l/r ) ( 8rT x /3r) are of approximately the same magnitude. Thus, the flow 
becomes more heavily influenced by the pressure gradient than by the turbu- 
lence gradients. In contrast with the situation for the upstream flow, the 

term 3T /3x in the downstream flow becomes more significant. Over most of 
xx 

the boundary layer, the magnitude of 3T xx /3x is about 20 percent of that for 
3p/3x, while it is about 25 to 30 percent of (l/r) ( 3rT rx /3r) in the lower half 
of the layer. 

The value of 3T xx /3x is still about 20 to 25 percent of 3p/3x at the 
downstream station (fig. 17(c)); while the shear stress gradient 
(l/r)(3rT rx /3r) tends to dominate nearly all of the boundary-layer flow. We 
see that the flow within and just downstream of the interaction is definitely 
not boundary- layer-like in that 3T /3x is not negligible; and, further, 
that the value of 3p/3x varies substantially across the boundary layer. 
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At the far downstream station (fig. 17(d)), the flow is more typical of 
that for a boundary layer, in that 3T^/ 3x has become quite small over most 
of the layer. 

Two additional points relative to the fluctuating flow can be made. In 
the Introduction, the susceptibility to separation when two interactions in 
the same boundary layer occur in proximity was noted. Two interactions were 
produced by a double cone similar to those used in reference 1 6 . The interac- 


tion resulting from the first shock wave produced boundary-layer separation. 

The second shock wave would also have been strong enough to produce separation 
of the undisturbed layer. However, when the second interaction occurred just 
downstream of the first (so that there were two distinct interactions), the 
second interaction did not produce a boundary-layer separation. This phenome- 
non could not be predicted by either of the methods of reference k or 5, but 
it may be explained with the aid of figure 1 6 as follows. It can be seen that 
the wall shear stress, T fx at the wall, continually increases throughout the 
region of adverse pressure gradient. This occurs because all of the turbulent 

mixing rates (e.g. , T and T ) have been increased so much as a result of 

xx rx 

the shock-wave - boundary-layer interaction that the momentum added to the flow- 
near the wall exceeds that removed in passing through the pressure rise. 
Therefore, if another shock-wave - boundary- layer interaction followed closely 
downstream of the one considered here, the boundary layer could sustain a 
larger pressure gradient than the initial layer without separating. 

Another point that is quite significant in the modeling of turbulent 


transport terms can be investigated with the aid of the shear stress T pre- 

rx 

sented in figure 1 6 . For so-called " equilibrium boundary layers,” a unique 
relationship exists between the shear stress and the mean velocity gradient 
3u/8y. The constant of proportionality between these two quantities is termed 
aii "eddy viscosity" or a "mixing length." We consider for simplicity here the 
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eddy viscosity model of the shear stress. The eddy viscosity e 


is defined 


as : 


One relationship proposed for 
(Clauser, ref. 52) is 


e = 


T /p 
rx 

3u/3y 


e for flow points away from the wall region 


e = const u 6* 
e 


which implies 


T /p = const u 6* — 
rx' K e 3y 


From the data in figure 1 6 and the "behavior of 6* (see fig. 10) and u g 
through the shock interaction region, we can assess the validity of this gen- 
eral type of modeling. We first note from table 1 that u & decreases by only 
about 10 percent between stations just upstream and just downstream of the 
interaction. From figure 10 we see that 6* decreases by about to percent. 

We may estimate that the average value of 3u/3y over the layer goes up by 
about 35 percent since u remains essentially unchanged and 6 decreases by 
about 35 percent. This would indicate that the quantity u g 6*(3u/3y) down- 
stream is about 70 percent of its value upstream. On the other hand, p 
increases by about a factor of 2 across the interaction. Therefore, if the 
equilibrium-shear-stress model were applicable, we would find T rx increasing 
by a factor of less than 2. We see from figure 1 6 , however, that T rx is 
nearly an order of magnitude larger over most of the boundary layer, even far 
downstream of the interaction (approximately nine downstream-boundary-layer 
thickness after the reflected shock system leaves the boundary layer). On the 
basis of this, we can conclude that boundary-layer flow downstream of the 
interaction is strongly out of equilibrium with the mean flow, a situation 
that should be accounted for by the shear-stress modeling used in analytic 
methods to describe this flow. Such a modeling has been proposed for use in 
compressible flows by Bradshaw and Ferris (ref. 53). 
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CONCLUSIONS 


The important conclusions derived from the present investigation are 
summarized below: 

1. The pressure rise to incipient separation is less than previously 
obtained in axially symmetric flow and is much less than corresponding planar 
two-dimensional flow. 

2. There are significant static-pressure and flow-angle variations across 
the boundary layer within and downstream of the interaction. 

3. The adverse pressure gradient significantly affects the total 
temperature distribution. 

U. The experimental turbulent stresses agree qualitatively with data 
obtained in incompressible flows with adverse pressure gradients. 

5. The Reynolds normal-stress term normally neglected in boundary-layer 
analyses is important within and just downstream of the interaction. 

6. The turbulent mixing rates are strongly out of equilibrium with the 
mean flow as a result of the interaction. 

7. Because of 6, novel modeling of the turbulent mixing rates is required, 
even far downstream of the interaction such as that proposed by Bradshaw and 
Ferris (ref. 53). 

Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif. 9^035, August 18, 1972 
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APPENDIX A 


EXPERIMENTAL APPARATUS 

The experimental investigation was conducted in the Heat-Power Laboratory 
of the Mechanical Engineering Building at the University of Washington, 

Seattle, Washington. The experimental apparatus used in this investigation 
consists of the supersonic wind tunnel (nozzle plus test section) and the 

instrumentation used to measure the flow within the test section. Figure 3 is 

a diagram of the test facility. 

The supersonic wind-tunnel apparatus is similar to that used by Seebaugh 
(ref. 15) , Mathews (ref. 4), and Teeter (ref. 1 6 ). However, significant 
improvements have been made in the ease of operation of the facility and in 
the actual aerodynamic components. A line sketch of the experimental facility 
is shown in figure 18. The wind tunnel is a continuous -flow facility using 
air supplied by two teflon-ring compressors capable of delivering in excess of 

1 lbm/s (2.2 kg/s) of dried and filtered air at 540° R (300° K) and TO psia 

( U . 82 x 10 5 N/m 3 ). The air passes through the first of two pressure regulators, 
through an electric heater to vary the air temperature, an additional dryer 
consisting of approximately 25 ft 3 (0.7 ni 3 ) of silica gel, a cartridge filter, 
the second pressure regulator, a quick- opening valve, a small gate valve, and 
into the tunnel plenum chamber. The pressure regulators allow continuous vari- 
ation of tunnel plenum (operating total) pressures from zero to the maximum 
pressure, depending on the size of tunnel in operation and the pressure level 
desired. The plenum chamber has a volume of approximately 3 ft 3 (0.08 m 3 ) and 
contains several screens and four 2-in. (5.07 cm) thicknesses of 1/4-in. 

(0.635 cm) cell honeycomb straightener , which ensures that the swirl is neg- 
ligible when the flow contracts from the plenum to the nozzle throat (area 
ratio approximately 150:1). 
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The plenum provides for the interchanging of the axially symmetric, 
plexiglass nozzle and test sections. The nozzle-test section facility that was 
used in the present study is shown in figure 19. The nozzle is 5.8 in. (lU.7 
cm) long and the constant-diameter test section is 6.5 in. (l6.5 cm) long. 

Also shown is the probe drive mechanism (the same one used by Grande, ref. ik) 
that drives either pitot probes or hot-wire anemometer probes in both the 
longitudinal (x) and radial (r) directions. The x and r locations can be 
set to the nearest 0.0005 in. (0.00127 cm) with the dial adjustments. The aft 
portion of the test section contains a 2. 2-in. -long (5.6 cm) by 0.2-in. -wide 
(0.51 cm) slot to allow for insertion and longitudinal translation of the 
probes. The contour of the nozzle is the same as that used by Seebaugh (ref. 
15) with the addition of a correction for the boundary-layer-displacement 
thickness. The coordinates of the nozzle are given in table 3. The downstream 
end of the test section attaches to the movable diffuser and model-support 
sting mechanism. An overall view of this test site, including plenum and 
diffuser, is shown in figure 20. 

The diffuser is supported by, and rigidly attached to, the carriage of a 
lathe bed. The carriage can be moved axially fore and aft to facilitate model 
changes in the tunnel and can be locked in position during operation. Also, 
the lathe bed is used to ensure precise alignment of the entire tunnel- 
diffuser combination each time the tunnel is opened and closed. The diffuser 
also contains the sting (see fig. 20 ) and sting translation mechanism. The 
sting is a case-hardened l/2-in.-diam. (l.27 cm) steel shaft running in three 
Thompson linear ball bushings. The axial motion of the sting is controlled 
externally by a pinion gear driving a rack gear on the sting. A counter indi- 
cates the sting location to within 0.005 in. (0.013 cm). 

The cone models used for generating the incident shock waves are mounted 
on the sting and may then be positioned anywhere on the centerline of the 
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tunnel. The cones themselves are straight cones turned from brass so that the 
base diameter is 1.30 in. (3.30 cm). The diffuser exhausts into a 150-ft 3 
O 4..25 m 3 ) vacuum tank maintained at a pressure of approximately 5 in. Hg 
(l.7 x 10 4 N/m 2 ) absolute, at maximum plenum pressure, by four Worthington 
Model 2 to steam ejectors. A view of the facility with the wind tunnel 
installed is shown in figure 21. 

A detailed survey of the flow within the nozzle and test section was con- 
ducted to define the mean- flow environment in which the boundary- layer inves- 
tigation took place and in which the hot-wire anemometer calibrations (appen- 
dix C) were made. This survey consisted of pitot-pressure measurements on the 
centerline from the nozzle throat through the test section, off-centerline 
pitot measurements in the test section, and surface-static pressures in the 
test section. 

The Mach number on the centerline deduced from these pitot pressure mea- 
surements is shown in figure 22. It can be seen that there is an overexpan- 
sion in the nozzle above the design value of 3.95. Recompression from this 
overexpansion then reduces the Mach number below the design value. Further 
expansions and compressions are evident. 

The surface-static-pressure distribution in the absence of a cone (fig. 
23) also indicates an overexpansion followed by a recompression. Near the 
origin of the test section, the wall pressure is low, and a definite rise 
occurs about 2 in. (5.1 cm) downstream. The off-centerline pitot measurements 
also confirm the presence of a recompression wave in this vicinity (about 
1 in. (2.5 cm) ahead of the region studied in the present investigation). 

Thus, the boundary layer on the nozzle-test section wall has a history of 
the favorable pressure gradient in the expansion section and of the slight 
adverse gradient induced by a recompression wave. 
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The off-centerline pitot measurements in the test section indicate that 
the Mach number external to the boundary layer has a value of 3.88 ± 0.02 over 
the length of the test section used in the present study. 

The static pressure taps used in this study were constructed as shown in 
figure 24(a), There are 33 taps located at 0.10-in. (0.254 cm) intervals in 
the main line of taps, plus nine symmetry taps located at 90° intervals around 
the periphery of the tunnel at longitudinal intervals of 1.00 in. (2.54 cm). 
These taps were used to confirm the symmetry of the flow when the conical shock 
generator was installed. 

The pitot pressures were measured with the probe shown in figure 24(b). 

The tip of the probe was made from 0.025-in. (0.063 cm) stainless steel tube 
with 0.006-in. (0.015 cm) walls. This tube was flattened and filed to produce 
a probe tip approximately 0.010 in. (0.025 cm) high with a 0.004-in. (0.010 cm) 
opening. The remainder of the probe is constructed of telescoping stainless 
tube as shown. 

The orifice dam used in the detection of incipient separation (appendix 
B) was located 1 in. (2.54 cm) downstream of one of the last symmetry taps. 

The dam itself was constructed as shown in figure 25. The dam is 0.010 in. 
(0.025 cm) wide and 0.005 in. (0.013 cm) high and extends along the tunnel 
periphery for 1/2 in. (1.27 cm). These dimensions give a length-to-height 
ratio of 100:1 so that the flow near the center of the dam where the pressure 
taps are located can be expected to be reasonably two-dimensional. 

The static pressures and the orifice dam pressures were measured with a 
CEC Model 4-312-0002, 0- to 10-psi absolute-pressure transducer. The absolute 
zero calibration point was obtained by a vacuum pump holding a pressure of 
approximately 60 to 80 ym Eg. The pitot pressures were measured with a CEC 
Model 4-312-0001, 25-psi differential-pressure transducer, with atmospheric 
pressure as reference. A Model J9D Scanivalve was used to aid in reading the 
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static pressures. All millivolt readings from the pressure transducers were 
obtained on a Dynasciences Model 440 digital multimeter whose lowest digit is 
1 yV. 

The instrumentation used for the fluctuating flow measurements consisted 
of the hot-wire probes and the anemometer system itself. A normal wire probe, 
used as a resistance thermometer, was used to measure total temperature. 

The hot-wire probes were constructed as shown in figure 26. The wire 
itself was of 5-ym diam. , etched tungsten wire supplied by the Sigmund Cohn 
Corporation, Mt. Vernon, New York. The wire supports (prongs) were made from 
high carbon steel, five-sided jeweler’s broaches. They were tapered from 
approximately 50 ym at the tip to approximately 200 ym 1/h in. (0.63 cm) from 
the tip. These prongs were set in a ceramic twin- conductor holder with epoxy 
resin so that approximately 1/8 in. (0.32 cm) of the prongs were exposed. The 
holder was approximately 0.0^0 in. (0.10 cm) in diameter and was wedge-shaped 
to minimize the influence of the probe support. Small twin-lead wire was used 
to connect the probes to the anemometer’s coaxial cable. The resulting probes 
had a distance between the probe tips of approximately 500 ym, and, thus, an 
i/d = 100. This is a lower value of i/d than is generally used and is the 

reason that individual calibration of each wire was required (see appendix C). 

The length was chosen to be small in order to obtain data as close to the wall 

as possible when the probe is in the r-x plane. With such a small distance 

between the probe tips, it was not feasible to use copper-plated wire ends 
(resulting in an even smaller sensing length of wire) similar, for example, to 
those used by Grande (ref. lb). Since tungsten cannot be soldered, the wires 
had to be arc welded to the prongs. During the present study, a micro-arc 
welder (fig. 27) was constructed and considerable time was spent learning how 
to attach wires without destroying the wire (and, at times, the entire prong!). 
The technique that was finally successful was to place the welding rod (a 
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0.003-in. (0.008 cm) silver wire) on top of the wire and in contact with the 
prong. The weld button was then depressed, causing a fine (roughly 10 ym) 
layer of silver to be deposited on the prong. Microscope examination of the 
weld showed that the wire was firmly embedded in the silver. 

The sketches of the probes shown in figure 2 6 indicate that the wires 
were mounted slack. This was done to prevent inordinate amounts of strain- 
gaging in the wire. No quantitative measure of the slackness can be given 
except to say that the two wires used for the final results of the present 
study contained strain-gage signals that were always less than 2 percent of 
the rms voltage. The only way found in the present study to ensure that the 
wires did not produce too much strain-gage signal was simply to construct and 
try them. Wires that appeared visually identical could give entirely differ- 
ent extraneous signals. Although the wires were mounted in a somewhat arbi- 
trary manner, the individual calibration of each wire allowed for any possible 
adverse effect of wire slack. 

The probes constructed as above, one a normal wire and the other yawed, 
lasted the entire test. The normal wire had a nominal cold (5^0° R (300° K) ) 
resistance of 1.03 Q and was run for 30 hr in the supersonic stream. The 
yawed wire had a nominal cold resistance of 1.26 Q and was run for U6 hr in 
the supersonic stream. During these run times the cold resistance changed 
only by about 3 percent (probably due to oxidation of the tungsten). Approx- 
imately 25 hot-wire probes were constructed and run during the preliminary 
phase of the present study; one should not get the impression that he need 
only construct two wires! 

The anemometer was a Disa model 55D01, constant-temperature system. The 
manufacturer’s quoted upper frequency response is 200 kHz, which corresponds 
to a decay time on the backside of a step input to the system of 0.80 ps . The 
decay-time oscilloscope traces for the wires constructed in the present study 


are shown in figure 28. Figures 28(a) and (b) show" decay traces for the 
normal wire at the minimum and maximum mass flows encountered in the study. 
Figure 28(c) shows the decay trace for the yawed wire at minimum mass flow. 

It can be seen that the decay time (the time required for the signal to 
decrease to within l/e of its final value) was always less than or equal to 
1 ys, indicating that the frequency response was in excess of 200 kHz through- 
out the range of flow. The blip in these time traces at about 3 to U ns indi- 
cates a slight overcompensation of the anemometer. The traces are shown for 
an overheat ratio a^ = 0.8. The frequency response decreases with decreasing 
overheat ratio; however, for the range of overheats used in the present study 
the frequency response was at least 200 kHz. 

This high value of frequency response permits resolution of disturbances 
in the boundary layer whose streamwise dimensions are greater than 6/5 in the 
outer part of the layer and greater than 6/8 in the lower part of the layer. 
Consider a "chunk” of self-correlated fluid moving in the boundary layer at a 

velocity u . For a turbulent boundary layer, u is approximately 0.6u 
c c e 

(Kistler and Chen, ref. 3*0, while in the inner part, the observed convection 
speeds may be as low as 0.4u e ( Vrebalovich , ref. 35). We may then use the 
equation relating u c and the upper frequency response f Q to find the mini- 
mum size of disturbance resolved by a given system: 

X . = u /f 

mm c o 

Taking f to be 200 kHz and u - 2200 ft/sec (670 m/s), X . - 0.0066 ft 

o e min 

(0.20 cm) in the outer part of the layer, while in the lower part 

X . =* 0.0 Ohb ft (0.13 cm), 

min 

To resolve a "chunk” of fluid passing the wire we need only consider X/2 
from the above equation. Thus, we can resolve disturbances larger than 0.0033 
ft (0.10 cm) in the outer part of the layer and 0.0022 ft ( 0.06 cm) in the 
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inner part. These values correspond, respectively, to 5/5 and 6/8. Kistler 
(ref. 29) estimates that 90 percent of the turbulent energy is contained in 
disturbances larger than 6 A, so that using the present anemometer and hot- 
wire combinations we should have been able to resolve over 90 percent of the 
turbulent energy. 
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APPENDIX B 


DETECTION OF INCIPIENT SEPARATION 

The strength of the shock wave used to generate the adverse pressure 
gradient imposed on the "boundary layer was near that required for incipient 
separation of the layer. This strength was chosen for two reasons. It is low 
enough that the large pressure fluctuations found in highly separated flows 
would not he present (such fluctuations significantly affect the interpreta- 
tion of the hot-wire anemometer signals). It is high enough to produce a 
pressure gradient that is realistic in terms of gradients that will exist in 
both internal and external flows for supersonic aircraft applications. 

Before a final decision was reached as to the strength of the shock (i.e., 
the cone generator angle), a brief experimental investigation was made of the 
incipient separation pressure rise for the flow of the present study. 

Seebaugh (ref. 15) concluded from his experimental results that the pressure 
rises for incipient separation in axially symmetric flow were much lower than 
those found by Kuehn (ref. 9, and substantiated by others) for planar two- 
dimensional flow. This appendix presents information on incipient separation 
pressure rise that confirms and/or extends the conclusions reached by Seebaugh 
(ref. 15) for axially symmetric flows. 

The flow conditions used in the present study were nearly the same as 
those used by Seebaugh (ref. 15). The Mach number was slightly higher, 3.88 
versus 3.78; and the boundary- layer- thickness Reynolds number Re^ slightly 
larger, 8.7 x 10 4 versus 6.2 x 10^. Both studies were conducted in nominally 
2-in. -di am. ( 5*08 cm) wind tunnels with the shock waves being generated by 
various angle cones placed in the center of the tunnel. 

Several methods for detecting the incipient separation condition have 
been proposed and used in the past. Six of these methods are listed below: 
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1. Determine whether a "knee" or "hump" occurs in the surface static 
pressure distribution. 

2. Determine whether the indicated reading of a pitot probe held on the 
surface becomes less than the surface static pressure in the absence of the 
pitot probe. 

3. Use measurements of pitot pressure near the wall and surface static 
pressure to find .Mach number and then extrapolate curves of Mach number versus 
y to zero Mach number (ref. l4 or 15). 

4. Use a light coating of oil on the surface and observe whether ridges 
of accumulation are present that could be associated with the stagnation 
points of separation and reattachment. 

5. Compare the static pressures observed just upstream and just down- 
stream of a small obstruction on the surface (the so-called orifice dam). 

6. Determine whether two reflected waves are present in the reflected 
shock system. 

Methods 2 and 3 are similar in that they neglect disturbance effects of the 
pitot probe as it is brought near the surface. A problem with this assumption 
was discussed in the Review of Previous Work. 

Methods 1, 2 (and essentially 3), 5, and 6 were investigated in the 
present study. Another method based on the introduction of alcohol at a 
minute rate into the boundary layer through a surface-static-pressure tap also 
was investigated. The physical principles involved with this technique are in 
some respects similar to the oil-flow technique (4); however, it is suffi- 
ciently different to consider it a new method. The alcohol injection method 
is discussed below. 

A tiny stream of alcohol is allowed to enter a surface-static-pressure 
tap at negligibly low velocity while the cone (and the shock-wave - boundary- 
layer interaction region) is positioned well ahead of the tap. The cone is 


then moved aft (see appendix C for details on motion of the conical generator) 
until the interaction region is quite close to the tap. On entering the tun- 
nel, the alcohol turns and flows downstream along the tunnel surface. Now, if 
the flow is separated, when the cone is moved just slightly further aft, the 
flow of alcohol will reverse and the alcohol will proceed entirely upstream. 
Without further motion of the cone, a well-developed ring of alcohol will 
propagate around the periphery of the tunnel wall. 

This change from downstream to upstream flowing alcohol is termed the 
"reattachment behavior," since the downs tream-to-upstre am change occurs at or 
very near what is most surely the reattachment point. So visually striking is 
the behavior that it is taken as the absolute criterion for the existence of 
separated flow. If no reattachment behavior can be observed, the boundary 
layer is assumed to have remained attached. The observed reattachment 
behavior is depicted schematically in figure 29. The location of the reattach- 
ment point can be determined quantitatively with a fair degree of certainty by 
the above procedure, since the presence of the alcohol itself probably alters 
the downstream to upstream behavior very little. 

The location of the separation point, however, cannot be quantitatively 
located with certainty since the ring of alcohol that forms in the upstream 
portion of the separated flow most likely alters the separation point from 
that which would exist with no alcohol ring present. The location of the 
upstream edge of the alcohol ring is probably farther forward than the separa- 
tion point in the absence of alcohol due to the buoyancy effect discussed in 
reference 10. 

Using the alcohol-injection technique, it was found that cones whose 
half-angles were 10° or larger separated the boundary layer, but that a 9° 
cone did not. In order to quantify this result somewhat, the results obtained 
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using the orifice dam (see appendix A) were compared with the results from the 
alcohol-injection study. 

The results of the orifice dam surveys made in the present study are 
shown in figure 30. The data were taken at x intervals of 0.015 in. (0.038 
cm) (Ax/ 6 - 0.07). The pressures from the upstream orifice and downstream 
orifice on the orifice dam are shown in comparison with the undisturbed static 
pressure (i.e., the surface static pressure in the absence of the dam). Also 
shown for the cases with separation are the locations of the reattachment 
point and the position of the leading edge of the alcohol ring. 

Figure 30 (a) shows the orifice dam results for a 12° half-angle cone. 

The reattachment point from the alcohol-injection technique matches exactly 
the point where the downstream pressure crosses the undisturbed pressure. 

There is some uncertainty in the location of reattachment due to uncertainty 
in the positioning of the cone relative to the pressure taps. An error band 
is indicated for the reattachment point. The position of the leading edge of 
the alcohol ring is also indicated. As previously noted, this is probably the 
farthest forward position at which separation could be occurring. Just down- 
stream of this location, the upstream and downstream orifice dam pressures 
cross. This indicates that the orifice dam is located in a region of total 
backflow - that is, a separated region. Thus the separation point must lie 
upstream of this crossover point. The leading edge of the alcohol is at the 
same location as the crossover between the downstream and undisturbed pres- 
sures. The separation point was chosen to be at the center of the extremes 
indicated above. This is shown in figure 30(a). The length of separation 

Z is indicated, 

sep 

Similar results are indicated for the 10° cone in figure 30 (b). For this 
case, a crossover of the upstream and downstream pressures does not occur. 
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The uncertainty in the separation point is assumed to be the same as in figure 
30(a). With this assumption, the value of & for the 10° cone is as shown. 

S 0JJ 

The results for the 9° cone are shown in figure 30(c). Reattachment 
behavior from the alcohol injection method was not observed for this case, so 
no alcohol locations are shown. The reattachment point (if any) was taken as 
the location of the downstream crossover point between the undisturbed and 
downstream pressures. The separation point was located by assuming that the 
forwardmost possible upstream location corresponds to the upstream crossover 
of the undisturbed and downstream pressures , and applying the uncertainty band 
as in figures 30 (a) and (b). A range of no separation to slight separation is 
indicated for the 9° cone. 

These separation lengths and their related uncertainties are shown in 
figure 31 as a function of cone half-angle. This plot indicates that a zero 
separation length (i.e. , the incipient separation case) could occur anywhere 
between 8.1° and 9.3°, the most likely value being 8.9°. This is in agreement 
with results from the alcohol injection, which showed that a 9° cone was the 
incipient separation case. Thus, flow with the 9° cone has been considered 
the incipient separation case in the present study. 

This value of incipient separation cone angle is lower than that found by 
Seebaugh (ref. 15) and corresponds to much lower pressure rises than Kuehn’s 
(ref. 9) planar values. These values are compared in figure 32. 

As noted earlier, other methods for determining incipient separation were 
investigated. Method 1, determining whether a hump occurs in the surface- 
static-pressure distribution, can be examined with the aid of figure 30 (a). 

For this case, a separated region existed but no hump in the undisturbed pres- 
sure distribution is evident. Thus, separation may occur without producing a 
hump in the pressure distribution. This was further confirmed in the 10° 
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interaction by taking undisturbed static-pressure data at every 0.0015 in. 
(0.0038 cm) , or one-tenth the spacing used in figure 30. 

Method 2 , and essentially method 3, can be tested by a comparison of the 
pitot pressure reading when the probe tip is at the wall with the undisturbed 
static pressure through an interaction known to be separated. The comparison 
for the 10° cone is shown in figure 33. Over the approximate separation 
length indicated, the pitot pressure might be expected to be less than the 
static pressure. This is not the case, however, so methods of detecting incip- 
ient separation based on pitot-pressure measurements near the wall are 
probably inaccurate. 

Finally, method 6, which bases the existence of separated flow on whether 
two reflected waves are present in the reflected shock system, was studied. 

As noted in the main body of the thesis, a well-developed, two-wave reflected 
shock system existed for the 9° cone interaction. For this interaction the 
flow was very near the condition for incipient separation. It is doubtful 
that this well-developed shock system could have suddenly appeared due to the 
presence of a slight (if any) separated region. Therefore, the two-wave 
shock-system criterion should not be used without some reservation. 
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APPENDIX C 


EXPERIMENTAL OPERATING PROCEDURES AND DATA REDUCTION 

The procedures followed in obtaining the mean- flow and fluctuating- flow 
data and in reducing these data are given in this appendix. Several prelimi- 
nary steps were taken before data were obtained. The diffuser was retracted 
and the desired conical shock generator installed on the sting. The diffuser 
was then locked in place and the downstream vacuum tank was evacuated. While 
the tunnel was under a vacuum, all the static-pressure instrumentation was 
checked to ensure against leaks. The pressure-transducer excitation voltage 
was set and the transducers zeroed to their respective reference pressures. 

The pressure regulators were set to give the desired operating total pressure 
(53.7 psia (3.7 * 10 5 N/m 2 ), nominally). The quick-opening valve to the plenum 
was then opened and supersonic flow established in the test section. 

The cone was then moved to a position that would cause the peak: pressure 
rise to occur near a set of symmetry pressure taps. The sting-bearing supports 
were then adjusted on the basis of symmetry of the wall static pressures to 
ensure that the cone was on the tunnel centerline and at zero incidence angle. 

The streamwise variations of the data obtained in the present study were 
obtained by two techniques. First, for the undisturbed boundary-layer studies, 
the cone was moved to the aft end of the tunnel. Data were obtained from the 
static pressure ports at locations 2.60, 2.80, and 3.00 in, (6.60, 7.11, and 
7.62 cm) downstream of the test section origin. The pitot pressures, total 
temperatures, and hot-wire data for the undisturbed boundary layer were also 
obtained at these stations. 

The second method for obtaining streamwise variation was to hold the 
instrumentation at a fixed x station and move the cone (and interaction 
region) past the instrumentation. This fixed instrumentation station was 
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located 3.00 in. (7.62 cm) downstream of the test section origin. Moving the 
cone in this manner allowed data to he obtained over a length of 2.2 in. (5.7 
cm) or approximately 11 boundary- layer thicknesses. An advantage of this tech- 
nique was that the profile measurements could be made 1 in. (2.5^ cm) ahead of 
the probe entry slot in the tunnel wall so that no interference from the slot 
was likely. When the cone was moved and the instrumentation held fixed, the 
flow entered the interaction at a slightly different value of Re^; however, 
this effect is quite small. 

A check on these two procedures was made to ensure that the data were 
essentially the same whether the cone was moved relative to one fixed static- 
pressure tap, or whether the cone was fixed and the various static-pressure 
taps were used to obtain the static-pressure distribution. The surface static 
pressures obtained by these two techniques agreed to within 2 percent every- 
where in the region in which data were obtained. 

All of the x stations referred to in the main body of the thesis were 
measured relative to the cone tip. A cone tip x station of 2.80 in. (7.11 
cm) corresponds to the tunnel x station of 3.00 in. (7.62 cm) (measured from 
the test section origin). The undisturbed boundary-layer data taken at tunnel 
x stations of 2.60 and 2.80 in. (6.60 and 7.11 cm) then correspond to cone 
tip x stations of 2.U0 and 2.60 in. (6.10 and 6.60 cm), respectively. Thus, 
data obtained at cone tip stations 2.U0, 2.60, and 2.80 in. (6.10, 6.60, and 
7.11 cm) provide information on the actual boundary- layer development ahead of 
the interaction, while data for cone-tip stations from x = 3.00 in. (7.62 cm) 
to x = U.60 in. (11.68 cm) show the effect of the interaction on the boundary 
layer holding the instrumentation station fixed and moving the cone. 

The correspondence between tunnel x station and cone-tip x station 
was obtained by sighting across the plastic tunnel along peripheral scribe 
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lines located at 1-in. (2.54 cm) intervals starting at the test section origin. 
The accuracy of this procedure was estimated to he ±0.002 in. (0.005 cm). 

As a means of stabilizing the operating total pressure and temperature, 

the wind tunnel was allowed to run for at least 30 min before data were taken. 

All electronic equipment, including the pressure transducers and the hot-wire 
anemometer units, was allowed to warm up at least 1 hr before data were taken. 

The pitot pressures were taken at each of the 12 stations indicated in 
figure 4. The probe tip was first brought to the wall, so that the midplane 
of the probe was 0.005 in. (0.013 cm) from the wall (half of the probe height). 
The probe was then moved 0.005 in. (0.013 cm) toward the center of the tunnel 
and the first data point obtained. Successive points were obtained at 0.010 
in. (0.025 cm) increments. The probe could be positioned in the y direction 
to within 0.0005 in. (0.0013 cm). The longitudinal position of the pitot 
probe was determined by sighting across the peripheral scribe lines to within 
an estimated accuracy of ±0.002 in. (0.005 cm). 

Total temperatures were obtained in the present study by using a normal 
wire, constructed exactly as the normal hot-wire probe shown in figure 26 (a), 
as a resistance thermometer. The wire probe was placed in a holder similar to 
that used to hold the pitot probe tip and was positioned with the same probe 
drive apparatus. The location of the wire in the flow was determined in a 
manner similar to that described for determining the pitot tip location, and to 
within about the same accuracy. 

The wire was calibrated in the wind tunnel with the cone retracted to its 
farthest aft position. The calibration was accomplished by using established 
recovery factor n versus Mach number curves for various Reynolds numbers and 
changing the total temperature of the stream with the electric heater (fig. l8) 
and changing the total pressure. The recovery curves obtained from the data 
summarized in reference 54 are shown in figure 34. The wire temperature could 
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then be calculated for each Reynolds number. At the same time, the wire 
resistance was measured. These measurements yielded a thermal coefficient of 
resistivity 04 that was valid for the probe in the presence of flow. The 
value obtained for the linear coefficient 04 in the expression 


= 1 + a x (T - T ) 

R . 1 w ref 

ref 

was 1.95 * 10 - 3 per °R (3.5 * 10 " 3 per °K) for T = 5^0° R (300° K). With 
this coefficient and measurements of the wire’s resistance at each point in 
the flow field, T^_ was calculated from the information in figure 3 ^. 

The fluctuating properties of the flow were obtained by a hot-wire ane- 
mometer. The hot-wire probes were installed on a probe shaft that was similar 
to that for the pitot probe. Positioning of the hot wires was accomplished by 
sighting across the peripheral scribe lines. The position of the yawed wire 
was taken to be at the center of the wire. 

The Disa 55D01 const ant -temperature anemometer unit was operated as sug- 
gested by the manufacturer. The input bias was adjusted to achieve a stable 
bridge operation over the range of overheat ratios used. The gain setting was 
set at position 9 for all cases. The high-frequency filter was set at position 
2. The low-frequency gain was set to the high position. The bridge ratio was 
20 : 1 . 


For each point in the flow, the wire's cold (recovery) resistance was 

read, the overheat ratio ( 0.1 ^ a < 1 . 0 ) selected, and the operating resis- 

w 

tance of the wire was set on the anemometer. Data were obtained at five over- 
heat ratios for each wire position used. At a given point in the flow, the 
normal wire was operated in only one position; however, the yawed wire was 
operated in four positions - two in the r-x plane and two in the x -0 plane. 
Fluctuating measurements were taken at a total of 137 flow- field points. 

Thus, just over 3^00 data points were taken. 
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For each of these 3^00 points, a square-wave frequency-response check was 
made by adjusting the inductance and capacitance controls in the anemometer to 
insure a high upper frequency response (see appendix A). 

After the frequency response check, the value of the mean bridge voltage 
E and the rms of the fluctuating voltage <E r > were recorded. The next 
operating resistance was then set and the process repeated. 

Each wire was individually calibrated before being used to obtain data. 
From appendix D we see that to obtain fluctuations of the flow properties from 
the voltage fluctuations for the normal wire, two sensitivities are required: 
the sensitivity to mass fluctuations Ae an< ^ sensitivity to total tem- 
perature fluctuations Aep^. For the yawed wire, the sensitivity to the v 
component Ae^ also is required. Direct calibrations of the required sensi- 
tivities were made over the range of Mach and Reynolds numbers encountered in 
the study. This range is shown in figure 35. In the calibration process, the 
wires were mounted on the sting on the centerline of the supersonic nozzle. 

The Mach number seen by the wire was varied by moving the probe into various 
centerline stations in the nozzle. The Reynolds number (and thus pu) was 

varied by changing the total pressure. To obtain the Ae sensitivity, a 

v 

holder that could position the wire at -10°, -5°, 0°, -+*5° * +10° incidence 
angles to the free-stream flow direction was mounted on the sting. 

The calibration of the normal wire was done as follows. The wire was 
positioned in the nozzle to give one of the desired calibration Mach numbers. 
For this position, the mean bridge voltage was read for various total pres- 
sures with constant total temperature. The results of these readings were 
then plotted as shown for two example Mach numbers in figure 36, The deriva- 
tive 3 in E/3 £n pu is needed to obtain Ae . From figure 36, 

pu 

3 £n E/3 in p^ may be obtained. At a fixed Mach number and total tempera- 
ture, 3 in p^ = 3 in pu. Thus, the desired term can be obtained directly 
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from plots similar to that of figure 36. The logarithmic relationship 
hetween E and p^ is not linear, so the slope must he determined as a function 
of Reynolds number. There was no significant effect of Mach number or overheat 
on 3 in E/3 jin pu. The values of 3 jin E/3 in pu for the normal wire as a 

function of the Reynolds number are shown in figure 37. 

An attempt was made to obtain Ae^ direct calibration of the quantity 

3 jin E/3 jin T^. This attempt was not successful, since when T^ was changed 
the total pressure also had to be changed to hold pu constant (as required). 
The total pressure simply could not be controlled accurately enough to ensure 
that the change in voltage was produced solely from a change in total 
temperature . 

The method used to calculate Aeij^ is based on the expression for Aep^ 

given in appendix D. The sensitivity Ae j ^ is made up of Ae plus two 

physical constants of the gas and one of tungsten, which are well known, plus 

a term involving 3 jin E/3 Jin R and R itself. The terms E and R , which 

w w w 

were found at the time of the actual data acquisition, are then curve fit as 

suggested by Grande (ref. lk) to give 3 in E/3 Jin R . Thus, all the informa- 

w 

tion used in calculating Ae^ was either obtained ahead of time, known from 

physical constants, or obtained in place. The Aeq*^ obtained in this manner 

is considered to be directly calibrated. 

The yawed-wire sensitivities to pu and T were calibrated in the same 

t 

manner as for the normal wire. The values obtained for 3 Jin E/3 Jin pu for 
the yawed wire were 0.65 of those for the normal wire over the calibrated 
range of flow variables. The value of 0.65 indicates that the effective yaw 
angle of the wire was approximately 40.5°. No attempt was made to measure the 
actual wire yaw angle, but it was intended to be nominally *+5°. 

The value of Ae^ for the yawed wire was also directly calibrated. The 
value of 3 Jin E/3<j> is required for Ae^. The value of <j> was varied by 
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varying the angle of incidence of the probe holder, which was mounted on the 
sting. The lowest reliable overheat ratio for the yawed wire was found to be 
0.2 (compared with Morkovin and Phinney's value of 0.15, ref. 26). The values 
of 3 in E/3<(> obtained from this calibration are shown in figure 38. There 
was no significant effect of overheat ratio on these values (above 0.2) and 
virtually no effect of Mach number. 

The data- re duct ion procedures used for the mean and fluctuating flow are 
discussed below. The mean-flow data consisted of pitot pressure and surface 
and flow- field static pressures, and the wire resistances used in determining 
the total temperature. The data- re duct ion program described in reference b 
was used in the input static pressure mode and was modified to calculate the 
total temperature from the input wire resistance. Tabulated values of n 
versus Re^. obtained from the curves of figure 3^ and the equation 



where is the wire resistance at T^, were incorporated into the program 

to calculate the total temperature. 

The fluctuating-flow data were reduced to usable form with computing pro- 
grams that embody the analyses given for both the normal and yawed wires in 
appendix D. The programs are listed on the following pages - the first is for 
the normal wire and the second is for the yawed wire. The yawed wire program 
was used for both the r-x plane and x-0 plane data, and w f takes the role 
of v ? when data for the x-0 plane are reduced. 

For examination of the derivatives of the fluctuating-flow data, the data 

were smoothed by least-squares-fitting parabolas through five adjacent points, 

which allowed for a smooth curve to approximate the data locally rather than 

globally. The derivatives were then obtained directly from the smoothed data 

by the properties of each local parabola. 
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The normal hot wire was also used to locate the shock-wave positions. 

This was done by observing the qualitative behavior of the anemometer output 
voltage on an oscilloscope. In the region of a shock wave or compression 
wave, the slight unsteadiness of the waves gave rise to a large fluctuating 
voltage. These large fluctuations were quite distinct from the boundary- layer 
turbulence in the outer half of the boundary layer (and, of course, in the 
inviscid flow). The locations of the incident and reflected shock systems 
were mapped out by this technique and were used to aid in the interpretation 
of the pitot-pressure profiles. 
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FORMAT STATEMENTS 


1001 FORMAT ( 8F 1 0 • 5) 
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I- 10 EBAR<1) MEAN BRIDGE VOLTAGE AT FIRST OVERHEAT RATIOpVOLTS 

II- 20 EBAR (2) »» SECOND '• 
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SLOCAL IS DLN E/DLN RHO U AT EACH RECAL 
SLOV IS DLN E/D PHI AT EACH RECAL 



RECALtSLOCALf AND SLOV ARE THE WIRE CALIBRATION DATA 
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ANGLE = ANGLE/57.29578 
EBARR = 0.0 
NWIRE = WIRENO 


NTEST * TEST 

NRUN = RUN 

RS = RB+RC 

00 20 1=1 » NOOH 

RW < I ) = AW(I)*RR*RR 






UJ 


UJ 







2 


X 







HH 


HH 







►- 









3 









O 


o 







x 


















a: 


UI 







UJ 


X 







o 


HH 







x 


x 







o 











X 







o 


o 







2 


X 







rvj 


X 







< 


UI 









o 







CO 


2 







< 


< 








z 

X 







a 

o 

o 







UJ 

X 








CO 

o 

V 







3 

» 

< 








X 

X 








UJ 








Z 

2 

o 







o 

H 

2 







HH 

OJ 

< 







K 

•* 








< 

vO 

Q 







-I 

H 

UJ 







O 

A 

K 







a. 

HH 

< 




CO 



a 

'— ■* 

X 




UJ 



UJ 

3 

GO 




H 



H* 

X 

HH 




K 



Z 

X 

-J 




HH 



HH 

o 

< 


HH 


> 




-J 

o 


w 


HH 



o 

CO 



3 


H* 



z 

#• 

CO 


X 


HH 



< 

1- 

M 


X 


CO 

> 



HH 



o 


2 

1- 


X 

u. 


o 

-J 


Ui 

HH 


3 

UJ 

o 

CO 

CO 


CO 

> 


* 

X 

-1 

>X> 

* 



1*4 


O 


UJ 

• 

• 


UJ 

H- 


o 

-1 

X 

o 

o 


x 

►-4 


-1 

< 


♦ 

o 


HH 

(O 



a 

o 

4-* 

fH 



2 


UJ 

o 

LO 

HH 

\ 



Ui 


_l 

-J 

yC 

w 



u. 

CO 



CO 

• 

3 

HH 


o 



< 


o 

X 






1- 

-J 


X 

X 


2 


o 


< 

UJ 

o 

< 


o 


• X 

< 

V 

3 

-1 

CD 


HH 


•“» O 


UJ 

3 

if) 

UJ 



3 

* O 

if) 

X 

< 




< 


H- 2 

HH 


> 

II 

II 

• 

-J 

1 

UI •> 


»— 




HH 

3 


X hh 

H- 

2 

UI 




O 

O 

II 

z 

HH 

X 

HH 

HH 

II 

-J 

X 

II hh 

HH 

< 

►- 


w 


< 

a: 


< 

H- 


3 

3 


o 


v- <— 4 

H 



X 

X 

HH 



Hi CO 


-I 


X 

UJ 

w 



u. 


-1 


o 

_J 

pH 



UJ o 


< 



UJ 

TC 



a o 


a 


CO 

3 

o 








*-H 

ru 








n 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


uauuu 


o u o a cj 

o o 



71 


CD 

* 

v£> 

♦ 

* 

* 


CO 

CL 

♦ 


>* 

GO 

£X 

CL 

< 


CL 31 

< CL 

CD 

UJ \ 

* — 

a *-* 

n ~ 


CD 

* 


CL 

I 











U_ 

*— 4 • 

to 











* -t 

• 










UJ 

X ♦ 

o 










> 

3 OJ 

♦ 










a 

CO * 











3 

* * 

1— I 










O 

vO ~ 












* »— 4 

3 










CO 

o ~ 

CL 










UJ 

H-I CL 

CL 










cc 

CO < 

O 










< 

UJ CD 

-J 










3 

CL UJ 

CO 










a 

* # 











CO 

»— 4 

* 


CO 









3 ^ 



o 








h- 

* CVJ 

*— 1 


• 








CO 

3 — 



o 








< 

CL 00 



* 

OJ — 



OJ 




UJ 

Z * 

< 

> 


* X 



# 




-J 

—1 • 

i 

h- 

M 

* CL 



* 


CL 



< OJ 

OJ 

•— i 









< 

•> w 

>3* — * 

> 


o 





O 



a v 


»-4 

X < 

CL O 





o 


CO 

CO x • 

• 

>- 

o ♦ 

CL -1 



QC 


-1 


*— 4 

UJ o —• 

o 3 

►H 

o S' 

< < 


3 

< 


< 



w o 

ac 

CO 

z o 

,-4 CD 


3 

CD 




-I 

JZ H 

— II CL 

z 

•> • 

♦ LU It 


► 

UJ 


It 

• 

o 

O * 

^0 o 

Ll) 

«— 4 *— 4 3d 

• 

04 


• 


<— ♦ 

CL 

Qu — - 

$ — * — 1 

CO 

It 

O II ^ 

1—4 

II 

II 

r-4 



O 

O II ^ 

* ►-* CO 


•-« II 

o -» 





•—4 

II 

CO 

CO 

^ * 



-T ~ ~ 

It 

1 


II 

w* 


-J 

-1 3 

~ k in 



3 


OJ 1— i 

*-4 


3 



n cl 

W s£> 


>t *-• 

II — Qd 



w 


CL 

vD 


J >r (L CL 0L 0- 

1 


C Z 

*— 4 

If 

o 

1— » 

Z 



-J o 

< Q • 

h- 

o *: 

“> to -J 


o 

CO 


-J 

f-4 


< O -J 2D — J O 


Q < 

3 UJ < 

3 

Q * 

UJ 3 

<t 

3 


O O CO UJ CO » 

<— 4 »— » 


<\J 


uuuu 


CJ o o u o 


uo uuu uuu 


72 


OELETT(I) = EBAR(I)/100.*SLOPTT(I) 


CALL TAINT (RECAL*SLOV*RET «SLOPV ( I ) »6*2»NER*EM0N) 
DELEV ( I ) = EBAR(I)/100.*SLOPV(I) 
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APPENDIX D 


MEASUREMENT OF TURBULENT TRANSPORT QUANTITIES IN A COMPRESSIBLE BOUNDARY LAYER 

A brief review of the state of signal interpretation from hot-wire anemom- 
eters was given in the body of the thesis. This appendix summarizes the impor- 
tant features of the interpretation techniques and extends the techniques for 
wire calibration and shear-stress measurement. The procedures for signal 
interpretation followed in the present study were strongly based on the works 
of Morkovin and Phinney (ref. 26), Morkovin (ref. 27), and Kovasznay (refs. 28 
and 30 ) . 

We begin by considering the heat loss from a heated wire held, possibly at 
a yawed angle, in a gas stream. We may write the energy balance for a correctly 
compensated wire (neglecting end losses) as: 

W = const I 2 R = Nu, (M,Re. ,6 ,<j>) • it ilk. • (T - n T ) = H (Dl) 

w w t t t w t 

where 0 = T /T, . Differentiating this expression logarithmically and using 
w t 

well-known relationships among Mach number, Reynolds number and velocity, 
density, and total temperature, and defining = 9 to k^/3 to T^ and 
m = 9 to y,/9 in T , we see that 

b U "C 

d in W = d in H (D2) 

where 

d to W = 2d to I + Kd to T = 2d in E + d to 

WWW 

and 


w 


(R + R ) 2 
_ w S _ 


(D3) 
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d to H = 


8 to Nu 

8 to Nu, 
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t ( 9 

n t " m t 8 to Re t “ 

9 to 0 0— n 

— 



+ 1 


1 

[ 1 8 to n 8 to n 

1 8 to Nu^ 

T wr 

2a 8 to M m t 8 to Re 

2a 9 to M 


^ J 
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d to T\ 


8 to Nu. n 8 to Nu. 

1 
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t ( 1 t 1 

1 9 to n j 9 to n 


9 to Re^_ a 9 to M t WT , 

|a 9 to M 9 to Re^ 



9 to Nu, 


1 9 to n 


9 Jin Re^ T vr 9 Jin Re. 

9 to Nu. 


d to p + 


9 Jin Nu 


t + e 


9 to 0 0-ri 


d to u 


d to T 

w 


1 9 to n 


L wr 


9(f) 


9 cf> 


d(j) 


(D4) 


where d(f) = v T /u. The + sign in the last term of equation (Dk) is associ- 
ated with a v* that increases the heat loss from the wire, as depicted in 
figure 39. 

If we invoke the constant temperature operation assumptions 
(dtoR =dtoT = 0 ) , we have 

2d to E = [ ]d to T + [ ]d to u + [ ]d to p + [ ] d<f) 

w U 



(D5) 


where the square brackets are those that appear in equation (DU). A substan- 
tial saving in time and effort during the data- acquis it ion phase can be made 
by noting that 

d to E = d Jin E 
w 

where E is the bridge voltage. This is true since E = E [ (R +R )/R ] and 

D w w s w 

d Jin R^ is assumed to be zero. Thus, we may divide equation (D5) by 2 and 


write 
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d E - T = 1 ) * tt) + W( 10 ° 1 t) * W( 10 ° tt) 

(d6) 

The fluctuation of the physical variables then is related to the fluctua- 
tion of the bridge voltage by 


A / «l \ / ir* 

— ) ± Ae 

P\ P/ vy u j 


E' = Ae^lOO =^) + AeJlOO \) + AeJlOO £=■) ± AeJlOO ^r] (D7) 


The terms Ae^ are the fluctuation sensitivities for the respective variables 
and, from equations (d 6) and (D7), they have the form 


Ae = 


E P NU t 1 


3 in n 


p 200 |j) in Re t i ¥r 3 in Re^j 


s Q to Nu, _ 3 £n Nu, 
Ae = + ~ t 


1 1 1 3 £n n 


u 200 (j) &n Re^ a 3 in M. x vr |a Hn M ' 3 £n Re_^ 


3 in n 


Ae Tt ~ 200 


3 to Nu 

n, + 1 - m, - — - — - — 

t t 3 in Re 


3 in Nu 


t + e 


+ _1 


t k. 

d in n 


3 in 0 0-n 


1 3 to Nu ‘ 


x wr |^2 a 3 £n M t 3 inReJ 2a 3 in M 


Ae v ~ 200 


\ ~ n 3 £n Nu “ 

1 3 in n t 


[Jwr 


3<() 


3cf) 


These fluctuation sensitivities are completely general in that they apply 
in either subsonic or supersonic flow for arbitrary Reynolds numbers and for 
yawed or normal (Ae^ = 0) wires. For supersonic Mach numbers, the complexity 
of these sensitivities is reduced considerably. It has been found (e.g., ref. 
26 or 5*0 that for M sin <f> > 1.2, all derivatives with respect to Mach number 
can be considered small, and neglected. With this condition, the sensitivities 
can be written as 
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Thus, for M sin 4> > 1.2, the relationship between the fluctuating voltage and 
fluctuations in the physical variables is just 


E' = Aem [100 =^) + Ae flOO + Ae (lOO ^ 

M T t/ pu \ pu / v \ 


We may evaluate the term { 3 in Nu^/3 in 0 + 9/ ( 0— n ) } in equation (D8b). 

Using Morkovin’s A = (l/2)(3 £n R /9 iln I ) and K = 3 in R /3 £n T , it can 
w w w w w 

be shown that { 3 to Efu,/3 to 9 + 0/(0-n)) = K[l + (l/A )]. Using the relation 

t w 

between I and E we see that 
w 



Now K is simply the logarithmic slope of R versus T^. For tungsten wire, 
the curve of R^ versus T^ was given by Grande (ref. lb ) . The logarithmic 
slope 3 £n E/3 in R^ is obtained from the data taken at the time of the 
actual fluctuating measurements. 

At this point we depart from the usual assumption (e.g. , ref. l4) that 


Re^ sin <j> > 20, which allows the variation of n with Re^ to be neglected so 

that, for example, the Ae sensitivity becomes 
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Ae 


pu 


E 

'3 Jin Nu ~ 

200 

3 Jin Re, 


L t-J 


This simplification would be desirable since, for very long wires, 

3 to Nu^/3 to Re^ = 1/2, and Ae pu thus becomes E/UOO. For act ual wires, 
however, the effects of end losses reduce this value, and approximations must 
be made to account for the influence of the losses. 

During the present study, it became apparent that the assumption 
Re^_ sin <{> > 20 would not be suitable for the boundary-layer measurements. To 
avoid the use of this assumption, therefore, each wire was directly calibrated 
over the entire range of Mach and Reynolds numbers encountered in the boundary 


layer. 

The direct calibration of the sensitivities is related to the above as 
follows. From equation (D9) we may write at once the mathematical form for 
the sensitivities 


A - -1- 9 E ''j 

T t 100 3 to T, 

X 

A = I 3 to E 

pu 100 3 to pu > 


(D10) 


A = E 8 £n E 
v " 100 

J 

The values of 8 £n E/8 in pu and 8 E/3^> were obtained by the calibra- 

tion method given in appendix C, where it also was noted that the value of 
8 in E/8 in T t could not be determined accurately. For Ae T , equation (D8b) 
was used with the substitutions 



in Nu 
x 

£n e 


+ 



2K 


8 £n E 

8 in R 
w 


+ 


R 

w 



and, following Morkovin, n = 0.885 and m = O.765. 

x t 

The resulting set of sensitivities therefore is obtained from terms that 
are either calibrated in place, calibrated beforehand, or are known physical 
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properties of air and tungsten. The only restriction on these sensitivities 
is that M sin 4> > 1.2, a condition satisfied in the present boundary- layer 
flow except for the first measuring station off the wall at two or three x 
locations. No end-loss corrections, peripheral calculations of Nusselt number, 
or similar calculations are required. We now discuss the use of these sensi- 
tivities in fluctuating flow measurements. We measure not E 1 but the rms of 
E f ; therefore, we must separate the (pu) 1 , Tj. , and v 1 contributions to the 
total rms E f . 

For a single normal wire in supersonic flow, Grande (ref. lb) has shown 
that the Morkovin diagram approach (ref. 27) and the Kovasznay diagram approach 
(ref. 28) are equivalent. The Kovasznay technique will be used here since it 
uses the basic sensing variables of the wire discussed above. This technique 
examines the mean square of E'. Thus, for normal wire Ae^ = 0, we have 


S 2 


e' 2 = r 2 (pyJ-L 2 - . 

(AeT t ) 2 (pu) 2 


(pu)’t; 

2r 1 

puT t 



10 4 


where r = -Ae /Ae^ . Now r changes with overheat ratio so that a plot of 

S versus r (the Kovasznay diagram) yields (pu) f 2 /(pu) 2 , (pu) 1 T^/puT^ , and 

T ,2 /T 2 as shown in reference 28. From these values and the definitions 
t t 


a 


1 

l+[(y-l)/2]M 2 


6 = (y-1)M 2 oi 


we can calculate 


> t 2 


(a+B) 2 



a 2 (pu) ' 2 + 2 
(a+B) 2 (pu) 2 


(pu)'T' 
a t 

(a+e ) 2 puT t 


u'T' 



(a+e ) 2 T 2 


a- 3 


(pu)'T' 


(a+e ) 2 puT t 


aB (pu) 1 2 
(a+e ) 2 (pu) 2 


and 
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p'u' = - — u'T' 

T 

For the output of a yawed wire, we examine the mean square of E' from 
the general equation for a yawed wire: 


10 0T' 

^2 = | Ae mm ur + __t + ioov 

mi mi -i-4- ^ r: 


FT = L 

pu pu 


L t T, 


< 4 %> 2 • 10,1 * ( Ae Tt > 2 \ • 104 - < 4 %> 2 — • 10 “ 

(pu ) 2 Tf u 2 

t 

(pu)'T' T’v f 

+2Ae Aerp • 10 4 ± 2Aern Ae • 10 4 

PU t PuT t t v T t u 

±2Ae Ae • 10 4 

pu v puu 


Conceptually, all the physical quantities required can he obtained from a 
single yawed wire operated at at least six overheat ratios. Attempts to do 
this in the present study with up to ten overheat ratios were unsuccessful 
because the equations were not sufficiently linearly independent. The 
following procedure was used instead. 

Mean square readings were taken before and after a 180° rotation about 
the probe axis. Then, taking the difference of these readings, we have: 


We let 


E 0° " E l80° " k Ae Tt Ae v “^T + Ae D u Ae 


v Ip - pu v puu 


(E 0° ~ E l 80 ° ] _ T t v ' _ r (pu) 1 v 1 

(l+Ae v Ae Tt • 10 4 ) T t u ^ uu 


Now a plot of S* versus r will immediately yield Tpu) 1 v T /pu 2 and T f v f /uT 

t t 

We can substitute these expressions, plus those for (pu)*/pu, T ! /T , and 

t t 
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(pu)’T' /puT obtained from the normal wire, back into equation (Dll) and 
"fc "b 

obtain v ,2 /u 2 . Further examination of the terms T|v' /T t u and (pu)'v'/pu 2 

leads to 

T t V * v'T' , „ u'v' 

— — = a + 6 — ~ 

T u uT u^ 

u 


and finally 


(pu) 1 v' _ u't' _ v'T 1 
pu 2 u 2 uT 


v’T 1 = uT 


T f v* 
1 t 


a+e 5 - 

t 


U ? V 1 = 


u 2 


T'v' 
1 t 

«+B f - 
t 


a+6 


a+6 


(pu) 1 v* 
pu 2 

(pu) ' v'~ 

pu 2 


We may also employ the yawed wire in the x-0 plane to obtain w' . 

The equations leading to w' 2 are identical to those leading to v , with 
the substitution w' = v' . Up to this point we have considered the mass-flux 
fluctuations and total-temperature fluctuations to be induced by vorticity 
and/or "temperature spottiness" and have neglected fluctuations in pressure (or 

sound) . 

When the hot wire is placed in a field in which only sound is present, the 
total-temperature and mass-flux fluctuation may be attributed to the sound 
field (ref. 28). The sensitivity of the wire to pressure fluctuations 
u = 100 (p' /yp) is: 


Ae^ = a(y-l) (l + n x M)Ae Tt + [l + ( n x / M ^ Ae pu 

where n x is the direction cosine of the wave. For example, if the sound wave 
is a Mach wave, then n x = -1/M, in which case, 

Ae = [1 - ( l/M 2 ) ] Ae 

7T P u 


8U 


Thus, for this case, it is related simply to (pu)'/pu by 

■„ = (P u )'/P^ 

1 - (l/M 2 ) 

Based on the findings of Grande (ref. l4), we conclude that the primary source 
of sound (aerodynamic noise) in the interaction is Mach waves. The measure- 
ments reported in the main body of the thesis were obtained using the above 
procedures. 
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APPENDIX E 


ACCURACY OF THE DATA 

The estimated accuracy of the mean- and fluctuating- flow data is discussed 
in this appendix. The quoted accuracies are "based on the past experience of 
others, the estimated uncertainty in the instrumentation of the present study, 
and, in some cases, the author’s opinion. The opinions are more in the nature 
of honest estimates rather than attempts to present glowing impressions of 
1 percent or 2 percent accuracy in one’s own work. The quoted accuracies "below 
are expressed as a percentage of the local value of each quantity. 

The estimated accuracies of the mean-flow data are for pitot pressure p^ , 
±1 percent; static pressure p, ±2 percent; and total temperature T , ±1/2 
percent. These uncertainties introduce other uncertainties in the derived 
mean- flow quantities such as u, p, T, etc, where the uncertainty is estimated 
to he less than 3 percent. The mean- flow measurements are quite "standard, 11 
and the 3 percent uncertainty in the individual terms should be considered 
extreme. Similarly, uncertainties exist in the products of these individual 
mean- flow terms - for example, pu and pu^ - that can be substantially larger 
than 3 percent, ±6 percent and ±9 percent, respectively. The last two uncer- 
tainties become significant if one attempts to check for a mass or momentum 
balance from the data. This is discussed further following a discussion of 
the uncertainty in the fluctuating- flow data. 

The estimated accuracies for the fluctuating- flow data are: <E’>, 

±2 percent; and E, ±1 percent. The errors in the reduced fluctuating data 
are, however, not simply reflections of these errors. The primary errors are 
introduced by the uncertainty in the sensitivity coefficients (appendix D) and 
the subsequent manipulations involving quantities of varying uncertainty. 
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The maximum errors in the pu- and v-sensitivity coefficients are esti- 
mated from Ae^^, ±5 percent from calibration; and Ae.^, ±5 percent from cali— 
brat ion* 

For the total temperature sensitivity, the following quantities, in addi- 
tion to the value of Ae^ u are required: K, ±3 percent (Grande, ref. l4); 

V ±1/2 P ercent ; V ±2 P ercent (Grande, ref. lb). These lead to an 
expected accuracy for the total temperature sensitivity Aeiji^_ of ±10 percent. 
This error is thought to he a maximum within the range of overheat ratios used 
in the present study , in spite of the fact that the errors at very low and very 
high overheat ratios can he much larger than 10 percent (Grande, ref. l4). 

With these errors, we can state the maximum expected uncertainty in the 
values of basic sensing variables (see appendix D) determined from the hot- 
wire measurements as listed below. These values should definitely be con- 
sidered the upper limits of the expected errors. 

<(pu) ? > , ±8 percent 

<T£> , ±12 percent 

(pu) f T£ , ±15 percent 

Tpu) 1 v ? , ±15 percent 

T£v f , ±15 percent 

<v*> , ±25 percent 

(pu) f w T , ±15 percent 

T£w f , ±15 percent 

<w'> , ±25 percent 

Examples of the extent to which these uncertainties may affect the 
reported values of Reynolds stresses are given below. Again these should be 
considered extreme. 
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T % ±15 percent 

xx 5 

T . ±20 percent 

rx 9 

T , ±30 percent 

rr 

, ±30 percent 

0 0 

T x9 , ±20 percent 

The derivatives of the Reynolds stresses reflect the uncertainties in the 
stresses themselves as modified hy the smoothing and differencing algorithm 
used to obtain the derivatives (appendix C). To estimate the uncertainty in 
the reported values of the derivatives, the following procedure was used. A 
10-percent variation was introduced in one value of T xx and in one value of 
T . The resulting variation of (l/r)(3rT /3r)and 3T /3x was then examined. 

rx 

A 10-percent variation of T rx introduced near a 10-percent variation in 
(l/r)(3rT rx /3r). A similar variation was found in ^ xx ^ x > so that the re P orted 
radial and axial derivatives carry uncertainties near those for the Reynolds 
stresses themselves. 

An attempt was made in the present study to verify the accuracy of the 
fluctuating flow measurements hy comparing the left- and right-hand sides of 
the integrated x-momentum equation. The x-momentum equation was integrated 
from the wall to the houndary-layer edge at each streamwise station. The con- 
tinuity equation was combined with the x-momentum equation to eliminate v. 
Both mean-flow terms on the left-hand side of the momentum equation were inte- 
grated numerically (Simpson’s rule) using the experimental data; terms on the 
right-hand side were treated in the same fashion. The streamwise derivatives 
that appear were determined hy the smoothing and differentiating procedure 
(appendix C) before integrating. The radial derivatives were integrated 
directly from the data so that no smoothing or differentiating was required. 
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This attempt was only qualitatively successful. Both terms on the left- 
hand side are of the order of 100. Their sum, which should be the right-hand 
side, is the order of 10 to 20, so that a 9-percent error in pu 2 as dis- 
cussed above could result in a value of this sum of as low as 0 and as high as 
20 to 30. On the right-hand side, the pressure gradient term is of the order 
of 5 to 10, the normal stress term is of the order of 1 to 3, and the shear 
stress term is of the order of 2 to 4. Thus, the right-hand side is of the 
order of 8 to IT* The left-hand side, being a small difference of two large 
numbers, simply cannot be known accurately enough to permit use of the inte- 
grated momentum equation as a suitable check on the accuracy of the data. 

The ratio of right- to left-hand sides, computed as discussed above, 
ranged randomly at the various streamwise stations from 0.7 to 2.5, indicating 
that the gradients of the experimental values of the Reynolds stresses are 
qualitatively correct; however, quantitative comparison cannot be realistically 
obt ained - 
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MEM FLOW DATA 
KEY FOR TABLE 1 
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PTINF 

TTINF 

PWALL 

PTP 
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x, cm 

y» cm 

p t ^ = 3.7 x 10 5 N/m 2 
Tf = 300° K 

u oo 

■^wall 

P tp 
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u, m/s 
p , kg/m 3 
T, °K 
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TABLE 2 


FLUCTUATING FLOW DATA 
KEY FOR TABLE 2 

X = x, cm 
Y = y, cm 

RHO'U' = ^T r u T , kg/m 2 -s 
RHO'V' = , kg/m 2 -s 

RHO'W' = p’w' , kg/m 2 -s 
RHOU'T' = kg °K/m 2 -s 

RHOV'T' = pv'T' , kg °K/m 2 -s 
RHOU ' U ’ = pu» 2 , N/m 2 
RHOV'V' = p7^, N/m 2 
RHOW'W' = pw' 2 , N/m 2 
RHOU ' V ' = pu'v' , N/m 2 
RHOU'W' = pu'w' , N/m 2 
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Reflected shock wave 



Figure 1.- Schematic of shock-wave - boundary-layer interaction region 








Figure k Location of measuring positions. 






































Reflected shock system 



Figure 6.- Experimentally deduced flow field. 





Incident shock wove Reflected shock system 



Figure 7.- Surface static pressures. 




















3. UO in. (8.6U cm), ^ wa ]_]_ 
Figure 8.- Continued. 





































































incident shock wave Reflected shock system 



Figure 12.- Boundary- layer mass flow 










Figure lU.- Turbulent mass transfer terms 





Shock position 



Figure lU.- Continued 






Figure lU.- Continued. 







cm 
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Figure lU # - Continued 



Figure lb m - Continued. 





Figure 1^.- Continued. 







Figure lU.- Continued. 








Figure lh.- Continued. 






Figure lU.- Continued. 







Figure 1*4.- Continued. 






Figure l4.- Continued. 



Figure lU.- Concluded 





Kg/m 0 -S 


(a) x = 2.U0 in. (6.10 cm) 

Figure 15*- Gradients of fluctuating terms from continuity equation 

(eq. (10)). 












Figure l6.- Reynolds stress terms. 



Figure 1 6 .- Continued. 





Figure 1 6 .- Continued 





Figure l6.- Continued 






Figure l6.- Continued 




Figure l6.- Continued. 






Figure l6.- Continued 





Figure l6.- Continued. 




Figure l6.- Continued 











Figure 1 6 .- Continued 





Figure l6,- Concluded, 





N/rn^ x 10 D 


(a) x = 2.U0 in. (6.10 cm) 

Figure IT.- Gradients of terms from right-hand side of x-momentum equation 

(eq. 11(a)). 
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Figure l8.- Line drawing of air flow through experimental apparatus. 
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Figure 20.- Experimental test site showing plenum, diffuser, and vacuum 



207 


Figure 21.- Experimental facility with nozzle/test section installed. 



'istribution of Mach number on nozzle /test section centerline. 





(a ) 

Note: All dimensions are in cm 


(a) Typical surface-static pressure instrumentation. 
Figure 2k.~ Pressure instrumentation. 



Note: All dimensions are in cm 


(l) ) Pitot-pressure probe. 
Figure 2U.- Concluded. 
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(a) Normal wire. 
Figure 2 6 .- Hot-wire probes. 
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Weld Welding 

button , electrode 



figure 27-- Micro-arc welder circuit. 
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Fl^drG 28 .- Hot — "W ire ciriSI!lOII10tS!I* UppSi' P^ggpi^Y'i oy 1 ypc^nriqo • a — P) 

50 mV/cm, 1 ys/cm. 


v 






















Figure 35-- Range of Mach and Reynolds numbers encountered in the present study 
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Figure 36.- Calibration technique for 3 to E/9 to pu 
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Figure 39.- Wire orientation for + sign convention in equation (j)h ) . 



